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Abstract 



£N| Type la supernovae (SNe la) play an important role in astrophysics and are crucial for the studies of stellar evolution, galaxy 
t-H evolution and cosmology. They are generally thought to be thermonuclear explosions of accreting carbon-oxygen white dwarfs 
(CO WDs) in close binaries, however, the nature of the mass donor star is still unclear. In this article, we review various progenitor 



(N 



models proposed in the past years and summarize many observational results that can be used to put constraints on the nature of 
their progenitors. We also discuss the origin of SN la diversity and the impacts of SN la progenitors on some fields. The cur- 



rently favourable progenitor model is the single-degenerate (SD) model, in which the WD accretes material from a non-degenerate 
companion star. This model may explain the similarities of most SNe la. It has long been argued that the double-degenerate 
OO (DD) model, which involves the merger of two CO WDs, may lead to an accretion-induced collapse rather than a thermonuclear 
explosion. However, recent observations of a few SNe la seem to support the DD model, and this model can produce normal SN 

i | la explosion under certain conditions. Additionally, the sub-luminous SNe la may be explained by the sub-Chandrasekhar mass 

rS { model. At present, it seems likely that more than one progenitor model, including some variants of the SD and DD models, may be 
required to explain the observed diversity of SNe la. 
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Type la supernova (SN la) explosions are among the most 
energetic events observed in the Universe. They are defined as 
those without hydrogen or helium lines in their spectra, but with 
strong Sill absorption lines around the maximum light (Filip- 
penko, 1997). They appear to be good cosmological distance 
indicators due to their high luminosities and remarkable unifor- 
mity, and thus are used for determining the cosmological pa- 
rameters (e.g. Qm and Q.\; Riess et al., 1998; Perlmutter et al., 
1999). This leads to the discovery of the accelerating expan- 
sion of the Universe that is driven by the mysterious dark en- 
ergy. SNe la are also a key part of our understanding of galactic 
chemical evolution owing to the main contribution of iron to 
their host galaxies (e.g. Greggio and Renzini, 1983; Matteucci 
and Greggio, 1986). In addition, they are accelerators of cosmic 
rays and as sources of kinetic energy in galaxy evolution pro- 
cesses (e.g. Helder et al., 2009; Powell et al., 201 1). The use of 
SNe la as standard candles is based on the assumption that all 
SNe la have similar progenitors and are highly homogeneous. 
However, several key issues related to the nature of their pro- 
genitors and explosion mechanism are still not well understood 
(Branch et al., 1995; Hillebrandt and Niemeyer, 2000). This 
may directly affects the reliability of the results of the current 
cosmological model and galactic chemical evolution model. 



Email addresses: wangboSynao .ac.cn (Bo Wang), 
zhanwenhanOynao . ac . en (Zhanwen Han) 



When SNe la are used as distance indicators, the Phillips 
relation is adopted, which is a phenomenological linear rela- 
tion between the absolute magnitude of SNe la and the magni- 
tude difference from its B-band maximum to 15 days after that 
(Phillips, 1993). The Phillips relation is based on the SN la 
sample of low redshift Universe (z < 0.05) and assumed to be 
valid at high redshift. This assumption is precarious since there 
is still no agreement on the nature of their progenitors. If the 
properties of SNe la evolve with the redshift, the results for 
cosmology might be different. In addition, more observational 
evidence indicates that not all SNe la obey the Phillips relation 
(e.g. Wang et al., 2006). 

Aside from the Phillips relation, many updated versions of 
this method are given to establish the relation between SN la 
intrinsic luminosities and the shape of their light curves. The 
stretch factor s method was proposed to measure the light curve 
shape by adjusting the scale on the time axis by a multiplica- 
tive factor (Perlmutter et al., 1997; Goldhaber et al., 2001). In 
addition, an empirical method based on multicolor light curve 
shapes has been developed to estimate the luminosity, distance, 
and total line-of-sight extinction of SNe la (Riess et al., 1998). 
Wang et al. (2005) presented a single post-maximum color pa- 
rameter AC12 (B — V color ~12 days after the B-band light max- 
imum), which empirically describes almost the full range of the 
observed SN la luminosities and gives tighter correlations with 
their luminosities, but the underpinning physics is still not un- 
derstood. Recently, Guy et al. (2005) used an innovative ap- 
proach to constrain the spectral energy distribution of SNe la, 
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parameterized continuously as a function of color and stretch 
factor s, and allow for the generation of light curve templates 
in arbitrary pass-bands. This method was known as the spec- 
tral adaptive light curve template method, which offers several 
practical advantages that make it easily applicable to high red- 
shift SNe la. The k-corrections are built into the model but not 
applied to the data, which allows one to propagate all the un- 
certainties directly from the measurement errors. 

It is widely accepted that SNe la arise from thermonuclear 
explosions of carbon-oxygen white dwarfs (CO WDs) in close 
binaries (Hoyle and Fowler, 1960; Nomoto et al., 1997). This 
hypothesis is supported by the fact that the amount of energy 
observed in SN explosions is equal to the amount that would 
be produced in the conversion of carbon and oxygen into iron 
(~10 51 erg; Thielemann et al., 2004). The energy released from 
the nuclear burning completely destroys the CO WD and pro- 
duces a large amount of 56 Ni. The optical/infrared light curves 
are powered by the radioactive decay of 56 Ni — » 56 Co — » 56 Fe. 
In order to trigger the carbon ignition, the mass of the CO WD 
must grow close to the Chandrasekhar (Ch) mass. When the 
WD increases its mass close to the Ch mass, it is thought to 
ignite near the center; at first the flame propagates subsonically 
as a deflagration, and in a second phase a detonation triggers, 
which propagates supersonically and completely destroys the 
CO WD (Hillebrandt and Niemeyer, 2000). The realistically 
conceivable way to make the WD grow to the Ch mass is via 
mass-transfer from a mass donor star in a close binary. How- 
ever, the nature of the mass donor star in the close binary is 
still uncertain, and no progenitor system before SN explosion 
has been conclusively identified. Additionally, there is some 
observational evidence that a subset of SNe la have progenitors 
with a mass exceeding or below the standard Ch mass limit (e.g. 
Howell et al., 2006; Foley et al., 2009; Wang et al., 2008a). 

Many progenitor models of SNe la have been proposed in 
the past years. The most popular progenitor models are single- 
degenerate (SD) and double-degenerate (DD) models. In Sect. 
2, we review various progenitor models, including some vari- 
ants of the SD and DD models proposed in the literature. We 
summarize some observational ways to test the current progen- 
itor models in Sect. 3, and introduce some objects that may be 
related to the progenitors and the surviving companion stars of 
SNe la in Sect. 4. We discuss the origin of SN la diversity 
and the impacts of SN la progenitors on some research fields 
in Sects. 5 and 6, respectively. Finally, a summary is given in 
Sect. 7. For more discussions on these subjects, see previous 
reviews on SN la progenitors (e.g. Branch et al., 1995; Hille- 
brandt and Niemeyer, 2000; Livio, 2000; Nomoto et al., 2003; 
Podsiadlowski, 2010; Maoz and Mannucci, 2012). 

2. Progenitor models 

2.1. Single-degenerate model 

In this model, a CO WD accretes hydrogen-rich or helium- 
rich material from a non-degenerate companion star, increases 
its mass to the Ch mass, and then explodes as a SN la (Whelan 
and Iben, 1973; Nomoto, 1982a). The SD model may explain 



the similarities of most SNe la, since SN la explosions in this 
model occur when the CO WD increases its mass to the maxi- 
mum stable mass (i.e. the Ch mass). In addition, the observed 
light curves and early time spectra of the majority of SNe la are 
in excellent agreement with the synthetic spectra of the SD Ch 
mass model (Nomoto et al., 1984; Hoflich et al., 1996; Nugent 
et al., 1997). 

The companion star in this model could be a main-sequence 
(MS) star or a subgiant star (WD + MS channel), or a red-giant 
star (WD + RG channel), or a helium star (WD + He star chan- 
nel) (Hachisu et al., 1996, 1999a,b; Li and van den Heuvel, 
1997; Langer et al., 2000; Han and Podsiadlowski, 2004; Fe- 
dorova et al., 2004; Meng et al, 2009; Wang et al., 2009a, 
2010a). The main problem for this class of models is that it 
is generally difficult to increase the mass of the WD by accre- 
tion. Whether the WD can grow in mass depends crucially on 
the mass-transfer rate and the evolution of the mass-transfer rate 
with time. (1) If the rate is too high, the system may enter into a 
common envelope (CE) phase; (2) if the rate is too low, the nu- 
clear burning is unstable that leads to nova explosions in which 
all the accreted matter is ejected. There is only a very narrow 
parameter range in which the WD can accrete H-rich or He- 
rich material and burn in a stable manner. This parameter range 
may be increased if the rotation affects the WD mass-accretion 
process (Yoon and Langer, 2004). 

An essential element in this model is the optically thick wind 
assumption, which enlarges the parameter space for produc- 
ing SNe la (Hachisu et al., 1996, 1999a,b; Li and van den 
Heuvel, 1997; Han and Podsiadlowski, 2004; Wang et al., 
2009a, 2010a). In this assumption, taking a MS donor star for 
an example, if the mass-transfer rate from the MS star exceeds 
a critical value, M cr , it is assumed that the accreted H burns 
steadily on the surface of the WD and that the H-rich material 
is converted into He at the rate of M cr . The unprocessed mat- 
ter is assumed to be lost from the binary system as an optically 
thick wind. However, this assumption is very sensitive to the 
Fe abundance, and it is likely that the wind does not work when 
the metallicity is lower than a certain valueQ 

2.1.1. WD + MS channel 

In the WD + MS channel (usually called the supersoft chan- 
nel), a CO WD in a binary system accretes H-rich material from 
a MS or a slightly evolved subgiant star. The accreted H-rich 
material is burned into He, and then the He is converted to car- 
bon and oxygen. When the CO WD increases its mass close 
to the Ch mass, it explodes as a SN la. Based on the evolu- 
tionary phase of the primordial primary (i.e. the massive star) 
at the beginning of the first Roche lobe overflow (RLOF), there 
are three evolutionary scenarios to form WD + MS systems and 



'At low enough metallicities (e.g. Z < 0.002), the optical depth of the wind 
would become small, and thus the wind-regulation mechanism would become 
ineffective (e.g. Kobayashi et al., 1998; Kobayashi and Nomoto, 2009). In this 
case, the binary system will pass through a CE phase before reaching the Ch 
mass. Thus, if this is true then there would be an obvious low-metallicity thresh- 
old for SNe la in comparison with SNe II. However, the metallicity threshold 
has not been found in observations (Prieto et al., 2008; Badenes et al., 2009a). 
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Figure 1 : Binary evolutionary scenarios of the WD + MS channel for producing SNe la. 



then produce SNe la (Fig. 1; for details see Wang et al., 2010a; 
also see Postnov and Yungelson, 2006; Meng et al, 2009). 

Scenario A: The primordial primary first fills its Roche lobe 
when it is in the Hertzsprung gap (HG) or first giant branch 
(FGB) stage (i.e. Case B mass-transfer defined by Kippenhahn 
and Weigert, 1967). In this case, due to a large mass-ratio or 
a convective envelope of the mass donor star, a CE may be 
formed (Paczyhski, 1976). After the CE ejection, the primary 
becomes a He star and continues to evolve. After the exhaus- 
tion of central He, the He star evolves to the RG stage. The He 
RG star that now contains a CO-core may fill its Roche lobe 
again due to the expansion itself, and transfer its remaining 
He-rich envelope onto the surface of the MS companion star, 
eventually leading to the formation of a CO WD + MS sys- 
tem. For this scenario, SN la explosions occur for the ranges 
M u ~ 4.0-7.0 M e , M 2>i ~ 1.0-2.0 M Q , and F ~ 5-30 days, 
where Mi i, M24 and P 1 are the initial masses of the primary 
and the secondary at zero age main-sequence (ZAMS), and the 
initial orbital period of the binary system. 

Scenario B: If the primordial primary is on the early asymp- 
totic giant branch (EAGB, i.e. He is exhausted in the center 
of the star while this star has a thick He-buming layer and the 
thermal pulses have not yet started), a CE may be formed due to 
the dynamically unstable mass-transfer. After the CE is ejected, 
a close He RG + MS binary may be produced; the binary or- 
bit decays in the process of the CE ejection and the primordial 
primary may evolve to a He RG that contains a CO-core. The 
He RG may fill its Roche lobe and start mass-transfer, which is 
likely stable and results in a CO WD + MS system. For this sce- 
nario, SN la explosions occur for the ranges M\ j ~ 2.5-6.5 M Q , 



M 2 ,i ~ 1.5 - 3.0 M and F ~ 200 - 900 days. 

Scenario C: The primordial primary fills its Roche lobe at 
the thermal pulsing asymptotic giant branch (TPAGB) stage. A 
CE is easily formed owing to the dynamically unstable mass- 
transfer during the RLOF. After the CE ejection, the primor- 
dial primary becomes a CO WD, then a CO WD + MS sys- 
tem is produced. For this scenario, SN la explosions occur 
for the ranges Mi; ~ 4.5-6.5 M Q , M 2 ,; ~ 1.5-3.5 M , and 
F > 1000 days. 

Among the three evolutionary scenarios above, models pre- 
dict that scenario A is the more significant route for producing 
SNe la (e.g. Wang et al., 2010a). The WD + MS channel has 
been identified in recent years as supersoft X-ray sources and 
recurrent novae (van den Heuvel et al., 1992; Rappaport et al., 
1994; Meng and Yang, 2010a). Many works have been con- 
centrated on this channel. Some authors studied the WD + MS 
channel with a simple analytical method to treat binary interac- 
tions (e.g. Hachisu et al., 1996, 1999a, 2008). Such analytic 
prescriptions may not describe some mass-transfer phases well 
enough, especially those occurring on a thermal time-scale. Li 
and van den Heuvel (1997) studied this channel from detailed 
binary evolution calculation with two WD masses (e.g. 1.0 
and 1.2 M ). Langer et al. (2000) investigated this channel 
for metallicities Z = 0.001 and 0.02, but they only studied 
Case A evolution (mass-transfer during the central H-burning 
phase). Han and Podsiadlowski (2004) carried out a detailed 
study of this channel including Case A and early Case B for 
Z = 0.02. The Galactic SN la birthrate from this study is 
0.6 - 1.1 x 10~ 3 yr~'. Following the studies of Han and Pod- 
siadlowski (2004), Meng et al. (2009) studied the WD + MS 
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channel comprehensively and systematically at various metal- 
licities. 

King et al. (2003) inferred that the mass-accretion rate on to 
the WD during dwarf nova outbursts can be sufficiently high to 
allow steady nuclear burning of the accreted matter and growth 
of the WD mass. Recently, Xu and Li (2009) also emphasized 
that, during the mass-transfer through the RLOF in the evolu- 
tion of WD binaries, the accreted material can form an accretion 
disc surrounding the WD, which may become thermally unsta- 
ble (at least during part of the mass-transfer lifetime), i.e. the 
mass-transfer rate is not equivalent to the mass-accretion rate 
onto the WD. By considering the effect of the thermal-viscous 
instability of accretion disk on the evolution of WD binaries, 
Wang et al. (2010a) recently enlarged the regions of the WD + 
MS channel for producing SNe la, and confirmed that WDs in 
this channel with an initial mass as low as 0.6 M can accrete 
efficiently and reach the Ch mass limit. Based on a detailed bi- 
nary population synthesis (BPS) approach]! they found that this 
channel is effective for producing SNe la (up to 1 .8 x 10 -3 yr _1 
in the Galaxy), which can account for about 2/3 of the observa- 
tions (see also Meng and Yang, 2010a). However, the parameter 
regions for producing SNe la in this model depend on many un- 
certain input parameters, in particular the duty cycle during the 
nova outbursts that is still poorly known. Additionally, whether 
dwarf nova outbursts can increase the mass of a WD is still a 
problem (e.g. Hachisu et al., 2010). 

2.1.2. WD + RG channel 

The mass donor star in this channel is a RG star, which is 
also called the symbiotic channel. There is one evolutionary 
scenario that can form WD + RG binaries and then produce 
SNe la (Fig. 2; for details see Wang et al, 2010a). Compared 
with the WD + MS channel, SNe la in the WD + RG channel 
are from wider primordial binaries. The primordial primary fills 
its Roche lobe at the TPAGB stage. A CE is easily formed 
due to the dynamically unstable mass-transfer during the RLOF. 
After the CE ejection, the primordial primary becomes a CO 
WD. The MS companion star continues to evolve until the RG 
stage, i.e. a CO WD + RG binary is formed. For the WD 
+ RG systems, SN la explosions occur for the ranges My ~ 
5.0-6.5 M , M 2ji ~ 1.0-1.5 M , andP 1 > 1500days. 

Unfortunately, the WD + RG binary usually undergoes a CE 
phase when the RG star overflows its Roche lobe. More im- 
portantly, the appropriate initial parameter space for producing 
SNe la in this channel is too small. Thus, WD + RG binaries 
seem to unlikely become a major way to form SNe la. Many au- 
thors claimed that the SN la birthrate via the WD + RG channel 
is much lower than that from the WD + MS channel (Yungel- 
son and Livio, 1998; Han and Podsiadlowski, 2004; Lii et al., 
2006; Wang et al., 2010a). The lowest initial WD mass in this 
channel for producing SNe la is about 1.0 M (e.g. Wang and 
Han, 2010a). In order to stabilize the mass-transfer process and 



2 BPS is a useful tool to simulate a large population of stars or binaries and 
can help understand processes that are difficult to observe directly or to model 
in detail (e.g. Han et al., 1995; Yungelson and Livio, 2000; Nelemans et al., 
2001). 
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Figure 2: Similar to Fig. 1, but for the WD + RG channel. 



avoid the formation of the CE, Hachisu et al. (1999b) assumed 
a mass-stripping model in which a stellar wind from the WD 
collides with the RG surface and strips some of the mass from 
the RG. They obtained a high SN la birthrate (-0.002 yr _1 ) for 
this channel. Here, Hachisu et al. (1999b) used equation (1) of 
Iben and Tutukov (1984) to estimate the birthrate, i.e. 

v = 0.2A<? AlogAyr 1 , (1) 

Jm a M z - 3 

where Aq, A log A, M\ and Mb are the appropriate ranges of 
the initial mass ratio, the initial separation, and the lower and 
upper limits of the primary mass for producing SNe la in units 
of solar masses, respectively. However, the birthrate is proba- 
bly overestimated, since some parameter spaces considered to 
produce SNe la in equation (1) may not contribute to SNe la. 

In symbiotic systems, WDs can accrete a fraction of the stel- 
lar wind from cool giants. It is generally believed that the stel- 
lar wind from a normal RG is expected to be largely spherical 
owing to the spherical stellar surface and isotropic radiation. 
However, the majority (>80%) of the observed planetary nebu- 
lae are found to have aspherical morphologies (Zuckerman and 
Aller, 1986). Additionally, the stellar winds from cool giants 
in symbiotic systems flow out in two ways: an equatorial disc 
and a spherical wind. In this context, by assuming an aspher- 
ical stellar wind with an equatorial disk from a RG, Lii et al. 
(2009) investigated the production of SNe la via the symbiotic 
channel. They estimated that the Galactic SN la birthrate via 
this channel is between 2.27 x 10~ 5 yr~' and 1.03 x 10~ 3 , and 
the theoretical SN la delay time (between the star formation and 
SN explosion) has a wide range from 0.07 to 5 Gyr. However, 
these results are greatly affected by the outflow velocity and the 
mass-loss rate of the equatorial disk. 
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The stellar wind from RG stars might be enhanced by tidal or 
other interactions with a companion. Tout and Eggleton (1988) 
brought the tidally enhanced stellar wind assumption to explain 
the mass inversion in RS CVn binaries. This assumption has 
been widely used to explain many phenomena related to giant 
star evolution in binaries (e.g. Han, 1998; van Winckel, 2003). 
The tidally enhanced stellar wind assumption has two advan- 
tages in the studies of symbiotic systems: (1) The WD may 
grow in mass substantially by accretion from stellar wind be- 
fore RLOF; (2) the mass-transfer may be stabilized because the 
mass ratio (M g i ant /M W D) can be much reduced at the onset of 
RLOF. By adopting the tidally enhanced stellar wind assump- 
tion, Chen et al. (201 1) recently argued that the parameter space 
of SN la progenitors can be extended to longer orbital peri- 
ods for the WD + RG channel (compared to the mass-stripping 
model of Hachisu et al., 1999b), and thus increase the birthrate 
up to 6.9 x 10~ 3 yr _1 , which is also probably overestimated due 
to the use of equation (1). Additionally, the parameter space of 
SN la progenitors strongly depends on the tidal wind enhance- 
ment parameter B w that is still poorly known. 

In a variant of the symbiotic channel, the mass-transfer from 
carbon-rich AGB stars with WD components can occur via stel- 
lar winds or RLOF (Iben and Tutukov, 1985). It has been sug- 
gested that an AGB donor star is in the progenitor system of SN 
2002ic, which is an atypical SN la with evidence for substantial 
amounts of hydrogen associated with the system (Hamuy et al., 
2003). Recently, Chiotellis et al. (2012) presented a WD with 
an AGB donor star for the SN remnant (SNR) of SN 1604, also 
known as Kepler's SNR. They argued that its main features can 
be explained by the model of a symbiotic binary consisting of 
a WD and an AGB donor star with an initial mass of 4-5 M . 
Detailed calculations of binary evolutionary model are needed 
to understand whether these WD components in WD + AGB 
binaries can result in SN la explosions. 

2.1.3. WD + He star channel 

A CO WD may also accrete helium-rich material from a He 
star or a He subgiant to increase its mass to the Ch mass, which 
is also known as the He star donor channel. There are three 
evolutionary scenarios to form WD + He star systems and then 
produce SNe la (see Fig. 3; for details see Wang et al., 2009b). 

Scenario A: The primordial primary first fills its Roche lobe 
when it is in the HG or FGB stage. At the end of the RLOF, 
the primary becomes a He star and continues to evolve. Af- 
ter the exhaustion of central He, the He star evolves to the RG 
stage. The He RG star that now contains a CO-core may fill its 
Roche lobe again due to the expansion of itself, and transfer its 
remaining He-rich envelope to the MS companion star, eventu- 
ally leading to the formation of a CO WD + MS system. After 
that, the MS companion star continues to evolve and fills its 
Roche lobe in the HG or FGB stage. A CE is possibly formed 
due to the dynamically unstable mass-transfer. If the CE can 
be ejected, a close CO WD + He star system is then produced. 
The CO WD + He star system continues to evolve, and the He 
star may fill its Roche lobe again (due to the orbit decay in- 
duced by the gravitational wave radiation or the expansion of 
the He star itself), and transfer some material onto the surface 



of the CO WD. The accreted He may be converted into carbon 
and oxygen via the He-shell burning, and the CO WD increases 
in mass and explodes as a SN la when its mass reaches the Ch 
mass. For this scenario, SN la explosions occur for the ranges 
M u ~5.0-8.0M o ,M 2 ,i ~ 2.0-6.5 M Q and F ~ 10-40days. 

Scenario B: If the primordial primary is on the EAGB stage 
at the onset of the RLOF, a CE may be formed due to the dy- 
namically unstable mass-transfer. After the CE is ejected, a 
close He RG + MS binary may be produced; the binary orbit 
decays in the procedure of the CE ejection and the primordial 
primary becomes a He RG. The He RG may fill its Roche lobe 
and start the mass-transfer, which is likely stable and results in 
a CO WD + MS system. The subsequent evolution of this sys- 
tem is similar to scenario A above, and may form a CO WD + 
He star system and finally produce a SN la. For this scenario, 
SN la explosions occur for the ranges My ~ 6.0 - 6.5 M e , 
M 2 ,i ~ 5.5 - 6.0 M and P { ~ 300 - 1000 days. 

Scenario C: The primordial primary fills its Roche lobe at 
the TPAGB stage, and the companion star evolves to the He- 
core burning stage. A double-core CE may be formed owing 
to the dynamically unstable mass-transfer during the RLOF. 
After the CE ejection, the primordial primary becomes a CO 
WD, and the companion star is a He star at the He-core burning 
stage, i.e. a CO WD + He star system is formed. The sub- 
sequent evolution of this system is similar to that in the above 
two scenarios, i.e. a SN la may be produced. For this scenario, 
SN la explosions occur for the ranges My ~ 5.5 - 6.5 M Q , 
M 2 ,i ~ 5.0 - 6.0 M and P i > 1000 days. 

SNe la from the He star donor channel can neatly avoid H 
lines, consistent with the defining spectral characteristic of most 
SNe la. Yoon and Langer (2003) followed the evolution of a 
WD + He star binary with a 1.0 M WD and a 1.6 M Q He star 
in a 0.124d orbit. In this binary, the WD accretes He from 
the He star and grows in mass to the Ch mass. Based on the 
optically thick wind assumption, Wang et al. (2009a) systemat- 
ically studied the He star donor channel. In the study, they car- 
ried out binary evolution calculations of this channel for about 
2600 close WD + He star binaries. The study showed the ini- 
tial parameter spaces for the progenitors of SNe la, and found 
that the minimum mass of CO WD for producing SNe la in this 
channel may be as low as 0.865 M . By using a detailed BPS 
approach, Wang et al. (2009b) found that the Galactic SN la 
birthrate from this channel is ~0.3 x 10~ 3 yr _1 and this chan- 
nel can produce SNe la with short delay times (~45-140Myr). 
Wang and Han (2010b) also studied the He star donor chan- 
nel with different metallicities. For a constant star-formation 
galaxy (like our own galaxy), they found that SN la birthrates 
increase with metallicity. If a single starburst is assumed (like 
in an elliptical galaxy), SNe la occur systematically earlier and 
the peak value of the birthrate is larger for a higher metallicity. 

2.2. Double-degenerate model 

In the DD model, SNe la arise from the merging of two close 
CO WDs that have a combined mass larger than or equal to 
the Ch mass (Tutukov and Yungelson, 1981; Iben and Tutukov, 
1984; Webbink, 1984). Both CO WDs are brought together by 
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Figure 3: Similar to Fig. 1, but for the WD + He channel. 



gravitational wave radiation on a timescale taw (Landau and 
Lifshitz, 1971), 



fcw(yr) = 8 x 10 7 (yr) x 



(M l + M 2 ) 
M\M 2 



1/3 



P 8/3 (h), 



(2) 



where P is the orbital period in hours, few in years and Mi, M 2 
in M Q . The delay time from the star formation to the occurrence 
of a SN la is equal to the sum of the timescale that the secondary 
star becomes a WD and the orbital decay time few- For the DD 
model, there are three binary evolutionary scenarios to form 
double CO WD systems and then produce SNe la, i.e. stable 
RLOF plus CE ejection scenario, CE ejection plus CE ejection 
scenario and exposed core plus CE ejection scenario (for details 
see Han, 1998). 

The DD model has the advantage that the theoretically pre- 
dicted merger rate is quite high, consistent with the observed 
SN la birthrate (e.g. Yungelson et al., 1994; Han, 1998; Nele- 
mans et al., 2001; Ruiter et al., 2009; Wang et al., 2OlOb)0 Im- 
portantly, this model can naturally explain the lack of H or He 
emission in the spectra of SNe la. As an additional argument in 
favor of the DD model, one may consider this model to explain 
some observed super-luminous SNe la (for more discussions 



3 Badenes and Maoz (2012) recently calculated the merger rate of binary 
WDs in the Galactic disk based on the observational data in the Sloan Digital 
Sky Survey. They claimed that there are not enough double WD systems with 
the super-Ch masses to reproduce the observed SN la birthrate in the context of 
the DD model. 



see Sect. 2.3.2). Furthermore, there are some double WD pro- 
genitor candidates that have been found in observations, and 
recent observations of a few SNe la seem to support the DD 
model (for more discussions see Sect. 3). However, the DD 
model has difficulties in explaining the similarities of most SNe 
la, since the merger mass in this model varies for different bi- 
naries and has a relatively wide range (~1.4 - 2.0 M ; Wang et 
al., 2010b). 

Most importantly, the merger of two WDs may result in an 
accretion-induced collapse to form a neutron star rather than 
a thermonuclear explosion (Nomoto and Iben, 1985; Saio and 
Nomoto, 1985; Timmes et al., 1994). In the process of a double- 
WD merger, once the less massive WD fills its Roche lobe, it 
is likely to be disrupted and rapidly accreted by the more mas- 
sive one. Meanwhile, the less massive WD is transformed into 
a disk-like structure around the more massive companion. It 
is usually assumed that in this configuration the temperature 
maximum is located at the "disk-dwarf" interface and that car- 
bon burning starts there. In this process, the carbon burning 
front propagates inward and then the CO WD is transformed 
into an O-Ne-Mg WD, which collapses to form a neutron star 
by electron capture on 24 Mg. 

There may be some parameter ranges where the accretion- 
induced collapse can be avoided (e.g. Piersanti et al., 2003; 
Yoon et al., 2007). Piersanti et al. (2003) suggested that the 
double WD merger process could be quite violent, and might 
lead to a SN la explosion under the right conditions. Pakmor 
et al. (2010) argued that the violent mergers of two equal - 
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mass CO WDs (~0.9 M , critical conditions for the success- 
ful initiation of a detonation) can be obtained, and may explain 
the formation of sub-luminous 1991bg-like objects. Although 
the light curve from the merger model is broader than that of 
SN 1991bg-like objects, the synthesized spectra, red color and 
low expansion velocities are all close to those observed for SN 
1991bg-like objects (Pakmor et al., 2010). In a further study, 
Pakmor et al. (201 1) claimed that a high mass-ratio is required 
for this model to work; for a primary mass of 0.9 M a mass- 
ratio of at least about 0.8. This result will affect the potential 
SN la birthrate of the DD model. We note that van Kerkwijk 
et al. (2010) came to a similar conclusion before Pakmor et 
al. (201 1), but that was in turn partially based on Pakmor et al. 

(2010) and Lorin-Aguilar et al. (2009). Adopting the results of 
Pakmor et al. (201 1) with a detailed BPS approach, Meng et al. 

(201 1) estimated that the sub-luminous events from this model 
may only account for not more than 1 % of all SNe la. 

Recently, by assuming that the moment at which the detona- 
tion forms is an artificial parameter, Pakmor et al. (2012) pre- 
sented a fully three-dimensional simulation of a violent merger 
of two CO WDs with masses of 0.9 M and 1.1 M Q , by combin- 
ing very high resolution and the exact initial conditions. They 
estimated that the simulation produces about 0.62 M of 56 Ni, 
and the synthetic multi-color light curves show good agreement 
with those observed for normal SNe la. Due to the small num- 
ber of such massive systems available, this model may only 
contribute a small fraction to the observed population of nor- 
mal SNe la. Future studies are needed to explore the parameter 
space of different WD masses and mass ratios in this scenario 
for normal SNe la, which is important in BPS studies. 

2.3. Potential progenitor models 

Besides the SD and DD models above, some variants of SD 
and DD models have been proposed to explain the observed di- 
versity of SNe la, such as the sub-Ch mass model, the super-Ch 
mass model, the single star model, the delayed dynamical insta- 
bility model, the spin-up/spin-down model, the core-degenerate 
model, the model of the collisions between two WDs, and the 
model of WDs near black holes, etc. 

2.3.1. Sub-Chandrasekhar mass model 

In this model, a CO WD accumulates a ~0.15 M e He layer 
with a total mass below the Ch mass (Nomoto 1982b; Woosley 
et al., 1986). In order to achieve the central densities necessary 
to produce iron-peak elements, the WD in this model needs a 
narrow mass range of ~0.9 - 1.1 M . The He may ignite off- 
center at the bottom of the He layer, resulting in an event known 
as Edge Lit Detonation (or Indirect Double Detonation). In 
this process, one detonation propagates outward through the He 
layer, while an inward propagating pressure wave compresses 
the CO core that ignites off-center, followed by an outward det- 
onation (e.g. Livne, 1990; Hoflich and Khokhlov, 1996). It 
is possible that sub-luminous 1991bg-like objects may be ex- 
plained by this model (Branch et al., 1995). Unfortunately, the 
sub-Ch mass model has difficulties in matching the observed 
light curves and spectroscopy of SNe la (Hoflich and Khokhlov, 



1996; Nugent et al., 1997), likely owing to the thickness of the 
He layer. 

Recently, Shen and Bildsten (2009) argued that, under some 
suitable conditions, a detonation in the WD might be achieved 
for even lower He layer masses than that in previous studies. 
By assuming that a detonation is successfully triggered in the 
He layer, Fink et al. (2010) claimed that the double detona- 
tions in sub-Ch mass WDs with low-mass He layers can be a 
robust explosion, leading to normal SN la brightness. Recent 
studies involving the sub-Ch mass WDs with subsequent nucle- 
osynthesis and radiative transfer calculations also indicate that 
the sub-Ch mass model could account for the range of the ob- 
served SN la brightness (Sim et al., 2010; Kromer et al., 2010). 
Additionally, BPS studies by Ruiter et al. (2009) predicted that 
there are a sufficient number of binaries with sub-Ch primary 
WDs to explain the observed birthrate of SNe la. However, it 
must be noted that it is difficult for the sub-Ch mass model to 
explain the similarities observed in most SNe la (e.g. Branch et 
al., 1995). 

2.3.2. Super-Chandrasekhar mass model 

The 56 Ni mass deduced from some SN la explosions strongly 
suggests the existence of super-Ch mass progenitors. SN 
2003fg was observed to be 2.2 times over-luminous than a nor- 
mal SN la, and the amount of 56 Ni was inferred to be 1.3 M Q , 
which requires a super-Ch mass WD explosion (~2. 1 M e ; How- 
ell et al., 2006). Following the discovery of SN 2003fg, 
three 2003fg-like events were also discovered, i.e. SN 2006gz 
(Hicken et al., 2007), SN 2007if (Scalzo et al., 2010; Yuan et 
al., 2010), and SN 2009dc (Yamanaka et al., 2009; Tanaka et 
al., 2010; Silverman et al., 2011). These super-luminous SNe 
la may raise the possibility that more than one progenitor model 
may lead to SNe la. 

It is usually assumed that these super-luminous SNe la are 
from the mergers of double WD systems, where the total mass 
of the DD systems is over the Ch mass. Meanwhile, a super-Ch 
WD may be also produced by a SD system, where the massive 
WD is supported by its rapid rotation, e.g. Maeda and Iwamoto 
(2009) claimed that the properties of SN 2003fg may be con- 
sistent with the aspherical explosion of a super-Ch WD, which 
is supported by its rapid rotation. Yoon and Langer (2004) ar- 
gued that WDs can rotate differentially for high mass-accretion 
rates of 3.0 x lO~ 7 M yr~'. By adopting the results of Yoon 
and Langer (2004), Chen and Li (2009) calculated the evolution 
of close binaries consisting of a CO WD and a MS star, and 
obtained the initial parameter space for super-Ch mass SN la 
progenitors. Within this parameter space, Meng et al. (2011) 
estimated that the upper limit of the contribution rate of these 
super-luminous SNe la to all SN la is less than 0.3%. Hachisu et 
al. (2012) recently made a comprehensive study of these super- 
luminous SNe la via the WD + MS channel, and suggested that 
these SNe la are born in a low metallicity environment as more 
massive initial CO WDs are required in this model. Meanwhile, 
Liu et al. (2010) also studied the He star donor channel to the 
formation of super-luminous SNe la by considering the effects 
of rapid differential rotation on the accreting WD. 
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Aside from the differential rotation, a super-Ch WD may also 
be supported by the WDs with strong magnetic fields due to 
the lifting effect. It has been found that ~10% of WDs have 
magnetic fields stronger than 1 MG (Liebert et al., 2003, 2005; 
Wickramasinghe and Ferrario, 2005). The mean mass of these 
magnetic WDs is ~0.93 M , compared with the mean mass 
of all WDs that is -0.56 M (e.g. Parthasarathy et al., 2007). 
Thus, the magnetic WDs are more easily to reach the Ch mass 
limit by accretion. The magnetic field may also affect some 
properties of SD progenitor systems, e.g. the mass-transfer rate, 
the critical mass-accretion rate and the thermonuclear reaction 
rate, etc. However, these effects are still unclear. Further studies 
are thus needed. 

2.3.3. Single star model 

Single star progenitor models have been considered by Iben 
and Renzini (1983) and Tout et al. (2008). In the absence of 
mass-loss, single massive star less than about 7 M will develop 
a degenerate CO-core when the star evolves to the AGB stage. 
The mass growing rate of the CO-core is controlled by the rate 
of the double shell burning. If the CO-core can grow to the 
Ch mass, it will produce a SN la. Under certain conditions, 
Tout et al. (2008) claimed that carbon can ignite at the center 
of the CO-core and the subsequent explosion would appear as 
a SN la. These single star progenitors are likely to be over 
2 M e , so this kind of SNe la should be associated with younger 
galaxies with recent star formation. The single star model was 
also proposed to explain the strongly circumstellar-interacting 
SN 2002ic (Hamuy et al., 2003). 

An important theoretical argument for this model is that the 
H-rich envelope in AGB star may be lost in a superwind be- 
fore the CO-core grows to the Ch mass, based on the envelope 
ejection criteria by Han et al. (1994) and Meng et al. (2008). 
Another problem for this model is that there should be far more 
SNe la than observed if a single star can naturally experience 
thermonuclear explosion. 

2.3.4. Delayed dynamical instability model 

This model is a variant of the WD + MS channel, which re- 
quires that the donor star is initially a relatively massive MS star 
(~3 M Q ) and that the system has experienced a delayed dynam- 
ical instability, resulting in a large amount of mass-loss from 
the system in the last a few 10 4 yr before SN explosion (Han 
and Podsiadlowski, 2006). The delayed dynamical instability 
model can reproduce the inferred H-rich circumstellar environ- 
ment, most likely with a disc-like geometry. Han and Pod- 
siadlowski (2006) claimed that the unusual properties of SN 
2002ic can be understood by the delayed dynamical instabil- 
ity model. Observationally, this model seems to be consistent 
with SN 2005gj (another 2002ic-like object) found by Nearby 
Supernova Factory observations (Aldering et al., 2006). 

However, in order for this model to be feasible, it requires a 
larger mass-accretion efficiency onto the WD than is assumed in 
present parametrizations. Based on a detailed BPS simulation, 
Han and Podsiadlowski (2006) estimated that not more than 1% 
SNe la should belong to this subclass of SNe la. Since this 
model requires an intermediate-mass secondary star, these SNe 



la should only be found in stellar populations with relatively 
recent star formation (e.g. with the last ~3 x 10 8 yr). 

2.3.5. Spin-up/spin-down model 

In the SD model, since the continued accretion of angular 
momentum can prevent the explosion of a WD, Justham (201 1) 
recently argued that it may be natural for the mass donor stars 
in the SD model to exhaust their envelopes and shrink rapidly 
before SN explosion, which may explain the lack of H or He 
in the spectra of SNe la, often seen as troublesome for the SD 
progenitor model. Di Stefano et al. (2011) also suggested that 
the CO WD is likely to achieve fast spin periods as the accreted 
mass carries angular momentum, which can increase the critical 
mass, M cr , needed for SN explosion. When the M a is higher 
than the maximum mass obtained by the WD, the WD must 
spin down before it explodes. This leads to a delay between 
the time at which the WD has completed its epoch of mass gain 
and the time of the SN explosion. However, the spin-down time 
is still unclear, which may have a large range from < 1 Myr to 
>1 Gyr (Lindblom, 1999; Yoon and Langer, 2005). The spin- 
down time may be important for the formation of the SNe la 
with long delay times. 

The spin-up/spin-down model may provide a route to explain 
the similarities and the diversity observed in SNe la. However, 
the birthrates, the delay times and the distributions of SN la 
explosion masses are still uncertain in this model. A detailed 
BPS studies are needed for this. 

2.3.6. Core-degenerate model 

Kashi and Soker (2011) recently investigated some possible 
outcomes of double WD mergers, in which these two compo- 
nents are made of CO. Most simulations and calculations of 
double WD mergers assume that a merger occurs a long time 
after the CE ejection, when these two WDs are already cold. 
In this model, Kashi and Soker (201 1) proposed that, a merger 
occurs within the final stages of the CE, whereas the CO-core 
is still hot. The merged hot core is supported by rotation until 
it slows down through the magnetic dipole radiation, and fi- 
nally explodes. Kashi and Soker (201 1) named this as the core- 
degenerate model, and claimed that this is another scenario to 
form a massive WD with super-Ch mass that might explode as 
a super-luminous SN la (see also Ilkov and Soker, 2012). A 
BPS study is required to determine the birthrate and delay time 
of this model, which are then compared with observations. 

2.3.7. Collisions of two WDs 

The WD number densities in globular clusters allow 
~10-100 times collisions between two WDs per year, and the 
observations of globular clusters in the nearby SO galaxy NGC 
7457 have detected a likely remnant of SNe la (Chomiuk et 
al., 2008). In this context, Raskin et al. (2009) explored colli- 
sions between two WDs as a way for producing SNe la. They 
carried out simulations of the collisions between two WDs 
(~O.6M ) at various impact parameters (the vertical separa- 
tion of the centers of the WDs). By taking impact parameters 
less than half of the WD radius before collision, they claimed 
that the SN explosions induced by such collisions can produce 
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~OAM Q of 56 Ni, making such objects potential candidates for 
sub-luminous SN la events. In a further study, Raskin et al. 
(2010) argued that two WD collisions could also realize super- 
Ch mass WD explosions (see also Rosswog et al., 2009a). 
However, this model predicts a very aspherical explosion, in- 
consistent with the small continuum polarization level in one 
of the observed super-luminous SNe la (i.e. SN 2009dc; see 
Tanaka et al., 2010). We note that collisions between two WDs 
are likely to happen in the dense environments of globular clus- 
ters, however the expected of which is still less frequent than 
that of the double WD mergers. 

2.3.8. WDs near black holes 

Wilson and Mathews (2004) proposed a new mechanism for 
producing SNe la, in which relativistic terms enhance the self- 
gravity of a CO WD when it passes near a black hole. They 
suggested that this relativistic compression can cause the cen- 
tral density of the WD to exceed the threshold forpycnonuclear 
reactions so that a thermonuclear runaway occurs. Dearborn 
et al. (2005) speculated that this mechanism might explain the 
observed 'mixed-morphology' of the Sgr A East SN remnant 
in the Galactic center. For more studies of this mechanism see 
Rosswog et al. (2008, 2009b). Due to the expected low rate of 
a WD passing near a black hole, the expected SN la birthrate 
from this mechanism should be significantly lower than that 
from normal SNe la. 

3. Observational constraints 

Many observational results can be used to constrain the SN 
la progenitor models, e.g. the properties of SN la host galaxies, 
the birthrates and delay times of SNe la, the candidate progen- 
itors of SNe la, the surviving companion stars of SNe la, the 
stripped mass of companions due to SN explosion, the signa- 
tures of gas outflows from some SN la progenitor systems, the 
wind-blown cavity in SN remnant, the early optical and UV 
emission of SNe la, the early radio and X-ray emission of SNe 
la, and the pre-explosion images and spectropolarimetry of SNe 
la, etc. 

3.1. SN la host galaxies 

There are some observational clues from the galaxies that 
host SNe la. SNe la have been known to occur both in young 
and old stellar populations (e.g. Branch and van den Bergh, 
1993), which implies that there is a time delay between the star 
formation and the SN explosion, ranging from much less than 
1 Gyr to at least several Gyr. In addition, SNe la in old popu- 
lation tend to be less luminous, and the most luminous SNe la 
appear to prefer young populations with recent star formation 
(Hamuy et al., 1996; Wang et al. 2008a). This indicates that the 
age of SNe la is an important parameter controlling at least part 
of SN la diversity. It was also established that super-luminous 
SNe la preferably occur in relatively metal poor environments 
with low-mass host galaxies, whereas sub-luminous SNe la oc- 
cur in non star-forming host galaxies with large stellar masses, 
such as elliptical galaxies (Neill et al., 2009; Taubenberger et 
al., 2011). 
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Figure 4: Evolution of Galactic SN la birthrates for a constant star formation 
rate (Z = 0.02, SFR = 5 M Q yr '). (From Wang et al., 2010b) 

The observational homogeneity of SNe la implies that a sin- 
gle progenitor system may produce most or all SNe la. How- 
ever, evidence for some observational diversity among SNe la, 
as well as evidence that SNe la can be produced by stellar pop- 
ulations that have a wide range of ages, raises the possibility 
that a variety of progenitor systems may be contributing. 

3.2. Birthrates of SNe la 

The observed SN la birthrate in our Galaxy is ~3xl0~ 3 yr -1 
(Cappellaro and Turatto, 1997), which can be used to constrain 
the progenitor models of SNe la. Based on a detailed BPS 
study, Wang et al. (2010b) systematically investigated Galac- 
tic SN la birthrates for the SD and DD models, where the SD 
model includes the WD + MS, WD + RG and WD + He star 
channels (see Fig. 4). They found that the Galactic SN la 
birthrate from the DD model is up to 2.9 x 10~ 3 yr _1 by as- 
suming that SNe la arise from the merging of two CO WDs 
that have a combined mass larger than or equal to the Ch mass, 
which is consistent with the birthrate inferred from observa- 
tions, whereas the total birthrates from the SD models can only 
account for about 2/3 of the observations, in which the birthrate 
from the WD + MS channel is up to 1.8 x 10~ 3 yr" 1 , the WD + 
RG channel is up to 3 x 10~ 5 yr 1 and the WD + He star channel 
is up to 0.3 x 10 3 yr -1 . The Galactic SN la birthrate from the 
WD + RG channel is too low to be compared with that of ob- 
servations, i.e. SNe la from this channel may be rare. However, 
further studies on this channel are necessary, since this channel 
may explain some SNe la with long delay times. In addition, it 
has been suggested that both recurrent novae, i.e. RS Oph and 
T CrB, are probable SN la progenitors and belong to the WD + 
RG channel (e.g. Belczyhski and Mikolajewska, 1998; Hachisu 
et al., 1999b; Sokoloski et al., 2006; Hachisu et al., 2007; Patat 
et al., 2007a, 2011). For other arguments in favour of the WD 
+ RG channel see Sects. 4.2 and 4.3. 

The SN la birthrate in galaxies is the convolution of the de- 
lay time distributions (DTDs) with the star formation history 
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(SFH): 

v(f)= f SFR(t-t')DTDs(t')dt', (3) 
Jo 

where SFR is the star formation rate, and t' is the delay time 
of a SN la. Due to a constant S FR adopted here, the SN la 
birthrate v(f) is only related to the DTDs, which can be ex- 
pressed by 

!0, t<h, 
DTDs'(t), h<t< t 2 , (4) 
0, t > t 2 , 

where t\ and t 2 are the minimum and maximum delay times of 
SNe la, respectively, and the DTDs' is the distribution of the 
delay times between t\ and t 2 . If t is larger than the t 2 , equation 
(3) can be written as 

rh 

v(f) = SFR DTDs'(t')dt' = constant. (5) 
Jt, 

Therefore, the SN la birthrates shown in Fig. 4 seem to be 
completely flat after the first rise. 

3.3. Delay time distributions 

The delay times of SNe la are defined as the time inter- 
val between the star formation and SN explosion. The vari- 
ous progenitor models of SNe la can be examined by compar- 
ing the delay time distributions (DTDs) expected from a pro- 
genitor model with that of observations. Many works involve 
the observational DTDs (e.g. Scannapieco and Bildsten, 2005; 
Mannucci et al., 2006, 2008; Forster et al., 2006; Aubourg et 
al., 2008; Botticella et al., 2008; Totani et al., 2008; Schawin- 
ski, 2009; Maoz et al., 2011). In recent years, three important 
observational results for SNe la have been proposed, i.e. the 
strong enhancement of the SN la birthrate in radio-loud early- 
type galaxies, the strong dependence of the SN la birthrate on 
the colors of the host galaxies, and the evolution of the SN la 
birthrate with redshift. Mannucci et al. (2006) claimed that 
these observational results can be best matched by a bimodal 
DTD, in which about half of SNe la explode soon after star- 
burst with a delay time less than lOOMyr, whereas others have 
a much wider distribution with a delay time ~3 Gyr. In a fur- 
ther study, Mannucci et al. (2008) suggested that 10% (weak 
bimodality) to 50% (strong bimodality) of all SNe la belong to 
the young SNe la. The existence of the young SN la popula- 
tion has also been confirmed by many other observations (e.g. 
Aubourg et al., 2008; Cooper et al., 2009; Thomson and Chary, 
2011), although with a wide range in defining the delay times 
of the young population. 

Maoz et al. (2011) presented a new method to recover 
the DTD, which can avoid some loss of information. In this 
method, the star formation history of every individual galaxy, 
or even every galaxy subunit, is convolved with a trial universal 
DTD, and the resulted current SN la birthrate is compared to the 
number of SNe la the galaxy hosted in their survey. They re- 
ported that a significant detection of both a prompt SN la com- 
ponent, that explodes within 420 Myr of star formation, and a 



delayed SNe la with population that explodes after 2.4 Gyr. Re- 
cently, a number of DTD measurements show that the DTD of 
SNe la follows the power-law form of r l (Maoz and Mannucci, 
2012). The power-law form is even different from the strong bi- 
modal DTD suggested by Mannucci et al. (2006), which might 
indicate that the two-component model is an insufficient de- 
scription for the observational data. We also note that there 
are many uncertainties in the observed DTDs, which are domi- 
nated by the uncertainties in galactic stellar populations and star 
formation histories (Maoz and Mannucci, 2012). 

Many BPS groups work on the theoretical DTDs of SNe la 
(e.g. Yungelson and Livio, 2000; Nelemans et al., 2001; Han 
and Podsiadlowski, 2004; Wang et al., 2009b, 2010a,b; Ruiter 
et al., 2009, 2011; Meng and Yang, 2010a; Mennekens et al., 
2010; Yu and Jeffery, 201 1; Claeys and Pols, 201 1). Other the- 
oretical DTDs of SNe la have been based on physically moti- 
vated mathematical parameterizations (e.g. Greggio and Ren- 
zini, 1983; Madau et al., 1998; Greggio, 2005, 2010). Re- 
cently, Nelemans et al. (201 1) collected data from different BPS 
groups and made a comparison. They found that the DTDs of 
different research groups for the DD model agree reasonably 
well, whereas the SD model have rather different results (see 
Fig. 5). One of the main differences in the results of the SD 
model is the mass-accretion efficiency with which the accreted 
H is added onto the surface of the WD (Nelemans et al., 201 1). 
However, the treatment of the mass-accretion efficiency cannot 
explain all the differences. Nelemans et al. are planning to do 
that in a forthcoming paper. For the SD model, Nelemans et al. 
(201 1) only considered systems with H-rich donor stars, not in- 
cluding the He-rich donor stars (Wang et al., 2009a). It is worth 
noting that the He star donor channel can produce SNe la effec- 
tively with short delay times (accounting for 14% of all SNe la 
in SD model; Wang et al., 2010b), which constitutes the weak 
bimodality as suggested by Mannucci et al. (2008). 

Hachisu et al. (2008) recently investigated new binary evo- 
lutionary models for SN la progenitors, with introducing the 
mass-stripping effect on a massive MS companion star by winds 
from a mass-accreting WD. This model can also provide a pos- 
sible way for producing young SNe la, but the model signifi- 
cantly depends on the efficiency of the artificial mass-stripping 
effect. Additionally, Chen and Li (2007) studied the WD + 
MS channel by considering a circumbinary disk which extracts 
the orbital angular momentum from the binary through tidal 
torques. This study also provides a possible way to produce 
SNe la with long delay times (~l-3 Gyr). 

3.4. Candidate progenitors 

3.4.1. Single- degenerate progenitors 

A number of WD binaries are known to be excellent candi- 
dates for SD progenitors of SNe la, e.g. U Sco, RS Oph and 
TCrB (Parthasarathy et al., 2007). All of these binaries contain 
WDs which are already close to the Ch mass, where the lat- 
ter two systems are symbiotic binaries containing a giant com- 
panion star (see Hachisu et al., 1999b). However, it is unclear 
whether these massive WD is a CO or an O-Ne-Mg WD; the 
latter is thought to collapse by forming a neutron star through 
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Figure 5: Rescaled delay time distributions from different BPS research groups, for the SD (left panel) and DD (right panel) models, respectively. The data points 
are taken from Maoz et al. (201 1). (From Nelemans et al., 201 1 ) 



electron capture on Mg rather than experience a thermonu- 
clear explosion (for more discussion see Sect. 4.2). Meanwhile, 
there are also two massive WD + He star systems (HD 49798 
with its WD companion and V445 Pup), which are good candi- 
dates of SN la progenitors. 

HD 49798 is a H depleted subdwarf 06 star and also a 
single-component spectroscopic binary with an orbital period 
of 1.548 d (Thackeray, 1970), which contains an X-ray pulsat- 
ing companion star (RX J0648.0-4418; Israel et al., 1997). The 
X-ray pulsating companion star is suggested to be a massive 
WD (Bisscheroux et al., 1997). Based on the pulse time delays 
and the inclination of the binary, constrained by the duration of 
the X-ray eclipse, Mereghetti et al. (2009) recently derived the 
masses of these two components. The corresponding masses 
are 1.50+0.05 M Q forHD 49798 and 1.28+0.05 M for the WD. 
According to a detailed binary evolution model with the opti- 
cally thick wind assumption, Wang and Han (2010c) found that 
the massive WD can increase its mass to the Ch mass after only 
a few 10 4 years. Thus, HD 49798 with its WD companion is a 
likely candidate of a SN la progenitor. 

V445 Pup is the first, and so far only, helium nova detected 
(Ashok and Banerjee, 2003; Kato and Hachisu, 2003). The 
outburst of V445 Pup was discovered on 30 December 2000 by 
Kanatsu (Kato et al., 2000). After that time, a dense dust shell 
was formed in the ejecta of the outburst, and the star became 
a strong infrared source, resulting in the star's fading below 
20 magnitudes in the V-band (Goranskij et al., 2010). From 
2003 to 2009, BVR observations by Goranskij et al. (2010) 
suggest that the dust absorption minimum finished in 2004, and 
the remnant reappeared at the level of 18.5 magnitudes in the 
V-band. Goranskij et al. (2010) reported that the most probable 
orbital period of the binary system is ~0.65 day. Based on the 
optically thick wind theory, Kato et al. (2008) presented a free- 
free emission dominated light curve model of V445 Pup. The 
light curve fitting in their study shows that the mass of the WD 
is more than 1 .35 M , and half of the accreted matter remains 
on the WD, leading to the mass increase of the WD. In addi- 



tion, the massive WD is a CO WD instead of an O-Ne-Mg WD, 
since no indication of neon was observed in the nebula-phase 
spectrum (Woudt and Steeghs, 2005). Therefore, V445 Pup is 
a strong candidate of a SN la progenitor (e.g. Kato et al., 2008; 
Woudt et al., 2009). 



3.4.2. Double-degenerate progenitors 

Several systematic searches for double WD systems have 
been made. The largest survey for this is SPY (ESO SN la Pro- 
genitor Survey; Napiwotzki et al., 2004; Nelemans et al., 2005; 
Geier et al., 2007), which aims at finding double WD systems 
as candidates of SN la progenitors. The only likely SN la pro- 
genitor in this sample is not a double WD system, but the WD 
+ sdB binary KPD 1930+2752 (Maxted et al., 2000). The or- 
bital period of this binary is 2.283 h, the mass of the sdB star 
is -0.55 M , and the mass of the WD is -0.97 M Q . The total 
mass (~1 .52 M ) and the merging time (<0.2 Gyr) of the binary 
indicate that it is a good candidate of a SN la progenitor (Geier 
et al., 2007). 

Recently, some other double WD systems have also been 
found, which may have the total mass close to the Ch mass, 
and possibly merge in the Hubble-time. These include a bi- 
nary WD 2020-425 with P olb ~ 0.3 day, M, + M 2 = 1.348 + 
0.045 M (Napiwotzki et al., 2007), V458 Vulpeculae with 
P 0lb ~ 0.068 day, M, ~ 0.6 M , M 2 > 1.0 M (Rodriguez- 
Gil et al., 2010), a close binary star SBS 1150+599A (double- 
degenerate nucleus of the planetary nebula TS 01) with P or b ~ 
0.163 day,M, = 0.54 + 0.02 M @ , M 2 ~ 0.86 M (Tovmassian et 
al., 2010), and GD687 that will evolve into a double WD sys- 
tem and merge to form a rare supermassive WD with the total 
mass at least 1 M (Geier et al., 2010). There are also some 
ongoing projects searching for double WD systems, e.g. the 
SWARMS survey by Badenes et al. (2009b) which is searching 
for compact WD binaries based on the spectroscopic catalog of 
the Sloan Digital Sky Survey. 
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3.5. Surviving companion stars 

A SN la explosion following the merger of two WDs will 
leave no compact remnant behind, whereas the companion star 
in the SD model will survive after a SN explosion and po- 
tentially be identifiable by virtue of its anomalous properties. 
Thus, one way to distinguish between the SD and DD models 
is to look at the center of a known SN la remnant to see whether 
any surviving companion star is present. A surviving compan- 
ion star in the SD model would evolve to a WD finally, and 
Hansen (2003) suggested that the SD model could potentially 
explain the properties of halo WDs (e.g. their space density and 
ages). Note that, there has been no conclusive proof yet that any 
individual object is the surviving companion star of a SN la. It 
will be a promising method to test SN la progenitor models by 
identifying their surviving companions. 

Han (2008) obtained many properties of the surviving 
companion stars of SNe la with intermediate delay times 
(100Myr-l Gyr) from the WD + MS channel, which are run- 
away stars moving away from the center of SN remnants. Wang 
and Han (2009) studied the properties of the companion stars of 
the SNe la with short delay times (<100Myr) from the He star 
donor channel, which are related to hypervelocity He stars (for 
more discussion see Sect. 4.5; also see Justham et al., 2009). 
Moreover, Wang and Han (2010d) recently obtained the proper- 
ties of the surviving companions of the SNe la with long delay 
times (>1 Gyr) from the WD + MS and WD + RG channels, 
providing a possible way to explain the formation of the popu- 
lation of single low-mass He WDs (<0.45 M ; for more discus- 
sion see Sect. 4.4; also see Justham et al., 2009). The properties 
of the surviving companion stars (e.g. the masses, the spatial 
velocities, the effective temperatures, the luminosities and the 
surface gravities, etc) can be verified by future observations. 

Tycho G was taken as the surviving companion of Tycho's 
SN by Ruiz-Lapuente et al. (2004). It has a space velocity of 
136 km/s, more than three times the mean velocity of the stars 
in the vicinity. Its surface gravity is log(g/cms~ 2 ) = 3.5 + 0.5, 
whereas the effective temperature is r e ff = 5750 ± 250K (Ruiz- 
Lapuente et al., 2004). These parameters are compatible with 
the properties of SN la surviving companions from the SD 
model (e.g. Han, 2008; Wang and Han, 2010d). However, 
Fuhrmann (2005) argued that Tycho G might be a Milky way 
thick-disk star that is coincidentally passing the vicinity of the 
remnant of Tycho's SN. Ihara et al. (2007) also argued that 
Tycho G may not be the companion star of Tycho's SN, since 
this star does not show any special properties in its spectrum; 
the surviving companions of SNe la would be contaminated by 
SN ejecta and show some special characteristics^! Recently, 
Gonzalez-Hernandez et al. (2009) presented some evidence 
that Tycho G may be enriched in 56 Ni, which could be the result 
of pollution of the atmosphere with the SN ejecta. 

By assuming that the companion star in the SD model is co- 
rotating with the binary orbit at the moment of the SN explo- 



Pan et al. (2012) studied the impact of SN la ejecta on MS, RG and He star 
companions with the FLASH code. They quantified the amount of contamina- 
tion on the companion star by the SN ejecta in their simulations, which might 
help to identify a companion star even a long time after the SN explosion. 



sion, the predicted rotational velocity of Tycho G is ~ 100 km/ s 
(e.g. Wang and Han, 2010d). However, the rapid rota- 
tion predicted by the SD model is not observed in Tycho G 
(7.5 + 2 km/s; Kerzendorf et al., 2009). This does not yet rule 
out that this star is the surviving companion. The inferred slow 
rotation of Tycho G may be related to the angular momentum 
loss induced by the rapid expansion of its outer shell. Recently, 
Pan et al. (2012) claimed that the post-impact companion star 
loses about half of its initial angular momentum for Tycho G, 
with the rotational velocity decreasing to a quarter of its initial 
rotational velocity, ~37 km/ s, which is closer to the observed 
rotational velocity (7.5 + 2 km/s). Therefore, whether Tycho 
G is the surviving companion of Tycho's SN is still quite de- 
batable. The confliction might be conquered by studying the 
interaction between the SN ejecta and the rotating companion 
star. 

We also note that Lu et al. (201 1) recently claimed that the 
angle between the direction of the non-thermal X-ray arc in Ty- 
cho's SNR to the explosive center and the proper motion veloc- 
ity of Tycho G is well consistent with the theoretical predictions 
and simulations. This supports Tycho G as the surviving com- 
panion of Tycho's SN. Lu et al. (201 1) also estimated the pa- 
rameters of the binary system before the SN explosion, which 
is useful for constraining progenitor models of SNe la. 

By investigating archival Hubble Space Telescope deep im- 
ages, Schaefer and Pagnotta (2012) recently reported that the 
central region of SNR 0509-67.5 (the site of a 1991T-like SN 
la explosion that occurred ~400 years ago) in the Large Mag- 
ellanic Cloud contains no surviving companion star. Thus, they 
argued that the progenitor of this particular SN la is a double 
WD system. In a subsequent work, Edwards et al. (2012) used 
the same method as in Schaefer and Pagnotta (2012) on SNR 
0519-69.0, which is a normal SN la remnant in the Large Mag- 
ellanic Cloud with an age of 600+200 years, and found that the 
99.73% error circle contains no post-MS stars for SNR 0519- 
69.0. Thus, Edwards et al. (2012) claimed to rule out the symbi- 
otic, recurrent nova, He star and spin-up/spin-down models for 
this particular SN. They argued that SNR 0519-69.0 might be 
formed from either a supersoft channel or a double WD merge. 
We note that, based on very short maximum spin-down times, 
Edwards et al. (2012) excluded the spin-up/spin-down model. 
However, if the spin-down time is much longer, the results in 
Edwards et al. (2012) might be different. 

3.6. Stripped mass of companions 

In the SD model, SN explosion will strip some mass of its 
non-degenerate companion star. By using two-dimensional Eu- 
lerian hydrodynamics simulations, Marietta et al. (2000) exam- 
ined the interaction of SN ejecta with a MS star, a subgiant star 
and a RG star. They claimed that the MS and subgiant com- 
panions lose ~ 10-20% of their mass after the SN explosion, 
and the RG companion loses about 96%-98% of its envelope. 
In this process, these stripped material is mixed with the SN 
ejecta. Since these stripped material is likely to be dominated 
by H, this should then lead to easily detectable H emission lines 
in the SN nebular phase. Unfortunately, no H has ever been de- 
tected in a normal SN la. The most recently observational up- 
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per limits on the amount of H detected are ~0.01 M (Leonard, 
2007)@ which may provide a strong constraint on the progen- 
itor model of SNe la. Additionally, based on the properties of 
the X-ray arc inside the Tycho's SNR, Lu et al. (201 1) also ob- 
tained a low stripped mass (<0.0083 M ), consistent with that 
from Leonard (2007). These observational limits are inconsis- 
tent with Marietta's predictions. 

Meng et al. (2007) used a simple analytical method to cal- 
culate the amount of the stripped masses. They obtained a 
lower limit of 0.035 M Q for the stripped mass, but their analytic 
method used oversimplified physics of the interaction between 
SN la ejecta and a companion star. Recently, many updated 
studies involve the effects of SN explosion on the companion 
star. However, more realistic stellar models for the compan- 
ion star do not show stripped mass as small as that close to the 
Leonard's observational limits, i.e. they do not resolve the con- 
flict between the theory and the observations (Pakmor et al., 
2008; Pan et al., 2010, 2012; Liu et al., 2012b). Thus, the high 
stripped mass from simulations may bring some problems for 
the SD model. The spin-up/spin-down model may explain the 
lack of H or He in SNe la (Justham, 2011; Di Stefano et al., 
2011; Hachisu et al., 2012). In addition, the mixture degree 
between the SN ejecta and the stripped material may also in- 
fluence the detection of H or He lines in the nebular spectra of 
SNe la. 

3. 7. Circumstellar material after SN explosion 

In the SD model, non-accreted material blown away from 
the binary system before SN explosion should remain as cir- 
cumstellar matter (CSM). Thus, the detection of CSM in SN la 
early spectra would support the SD model. Patat et al. (2007a) 
found some direct evidence on CSM in a normal SN la, i.e. SN 
2006X, which was also exceptional in its high ejecta velocity 
and high reddening (Wang et al., 2008b). Patat et al. (2007a) 
have observed a variation of Na I doublet lines immediately 
after the SN explosion, which is interpreted as arising from 
the ionization and subsequent recombination of Na in CSM. 
This strongly favours a SD progenitor for this SN. Patat et al. 
(2007a) suggested that the narrow lines may be explained by a 
recurrent nova. The time-variable Na I doublet absorption fea- 
tures are also found in SN 1999cl (Blondin et al., 2009) and 
SN 20071e (Simon et al., 2009). Patat et al. (2007a) argued 
that the CSM may be common in all SNe la, although there 
exists variation in its detect ability because of viewing angle ef- 
fects. However, in a subsequent work, Patat et al. (2007b) did 
not find the same spectral features in SN 2000cx as they did 
with SN 2006X, which indicates that there might be multiple 
SD progenitor models. Meanwhile, the derivation of smaller 
absorption ratio Ry (the ratio of the total to selective absorp- 
tion by dust) perhaps also suggests the presence of CSM dust 
around a subclass of SNe la (Wang et al., 2009c). 

More encouragingly, Sternberg et al. (2011) studied the ve- 
locity structure of absorbing material along the line of sight to 



'Leonard (2007) obtained deep spectroscopy in the late nebular phase of 
two well observed SNe la (SN 2005am and SN 2005cf), in search of the trace 
amounts of H and He that would be expected from the SD model. 



35 SNe la in nearby spiral galaxies via Na I doublet absorption 
features. They found a strong statistical preference for blue 
shifted structures, which are likely signatures of gas outflows 
from the SN la progenitor systems. They concluded that many 
SNe la in nearby spiral galaxies may originate in SD systems, 
and estimated that at least 20% of SNe la that occur in spiral 
galaxies are from the SD progenitors. Recently, Foley et al. 
(2012) reported that SNe la with blue shifted structures have 
higher ejecta velocities and redder colors at maximum bright- 
ness relative to the rest of the SN la population, which provides 
the link between the progenitor systems and properties of SN 
explosion. This result adds additional confirmation that some 
SNe la are produced from the SD model. However, Shen et al. 
(2012) argued that such gas outflow signatures could also be in- 
duced by winds and/or the mass ejected during the coalescence 
in the double WDs. 

3.8. SN remnants 

SN remnants (SNRs) are beautiful astronomical objects that 
are also of high scientific interest, since they provide direct in- 
sights into SN progenitor models and explosion mechanisms. 
Recent studies by Lu et al. (2011) suggested that the non- 
thermal X-ray arc in Tycho's SNR is a result of interaction be- 
tween the SN ejecta and the stripped mass of the companion, 
strengthening the motivation of studying the progenitor of a SN 
by studying its SNR. In addition, SNRs may reveal the metal- 
licity of SN progenitors (Badenes et al., 2008). 

Circumstellar matter (CSM) is predicted by the SD model, 
which was responsible for creating a low-density bubble (i.e. 
wind-blown cavity; Badenes et al., 2007). Its modification on 
larger scales will become apparent during the SNR phase. One 
of the obstacles the SD model faces is to search for this signa- 
tures from SNR observations. Badenes et al. (2007) searched 7 
young SN la remnants for the wind-blown cavities that would 
be expected in the SD model. Unfortunately, in every case it 
appears that the remnant is expanding into a constant density 
interstellar matter (i.e. there is no wind-blown cavity in these 
SN remnants). However, Williams et al. (2011) recently re- 
ported results from a multi-wavelength analysis of the Galactic 
SN remnant RCW 86 (remnant of SN 185 A.D.). From hydro- 
dynamic simulations, the observed characteristics of RCW 86 
are successfully reproduced by an off-center SN explosion in a 
low-density cavity carved by the progenitor system (Williams 
et al., 2011). This makes RCW 86 the first known case of a 
SN la in a wind-blown cavity. 

3.9. Early optical and UV emission of SNe la 

The presence of a non-degenerate companion in the SD 
model could leave an observable trace in the form of the optical 
and ultraviolet (UV) emission. Kasen (2010) showed that the 
collision of the SN ejecta with its companion should produce 
detectable optical and UV emission in the hours and days fol- 
lowing the SN explosion, which can be used to infer the radius 
of the companion. Thus, the early optical and UV observations 
of SN ejecta can directly test progenitor models. The optical 
and UV emission at early times forms mainly in the outer shells 
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of the SN ejecta, in which the unburned outer layers of the WD 
play an important role in shaping the appearance of the spec- 
trum. Kasen (2010) claimed that these emission would be ob- 
servable only under favorable viewing angles, and its intensity 
depends on the nature of the companion star. 

Hay den et al. (2010) looked for this signal in the rising por- 
tion of the B-band light curves of 108 SNe la from Sloan Digital 
Sky Survey, finding no strong evidence of a shock signature in 
the data. They constrained the companion in the SD model to 
be less than a 6 M e MS star, strongly disfavouring a RG star 
undergoing RLOF. Recently, Bianco et al. (2011) searched for 
the signature of a non-degenerate companion star in three years 
of SN Legacy Survey data by generating synthetic light curves 
accounting for the shock effects and comparing true and syn- 
thetic time series with Kolmogorov-Smirnov tests. Based on 
the constraining result that the shock effect is more prominent 
in rest-frame B than V band (for details see Fig. 3 of Kasen, 
2010), Bianco et al. (2011) excluded a contribution of WD + 
RG binaries to SN la explosions. However, a rather contradic- 
tory result for the shock effects was obtained by Ganeshalingam 
etal. (2011). 

These shock signatures predicted in Kasen (2010) are based 
on the assumption that the companion star fills its Roche lobe 
at the moment of a SN explosion. However, if the binary sep- 
aration is much larger than the radius of the companion star, 
the solid angle subtended by the companion would be much 
smaller. Thus, the shock effect would be lower. Justham (201 1) 
and Di Stefano et al. (2011) argued that the donor star in the 
SD model may shrink rapidly before the SN explosion, since 
it would exhaust its H-rich envelope during a long spin-down 
time of the rapidly rotating WD until the SN explosion. In 
this condition, the companion star would be a smaller target 
for the SN ejecta and produce a much smaller shock luminos- 
ity than the Roche lobe model considered in Kasen (2010) (see 
also Hachisu et al., 2012). Therefore, the early optical and UV 
emission of SN ejecta may be compatible with the SD model. 

In recent optical and UV observations, Wang et al. (2012) 
presented UV and optical photometry and early time spectra 
of four SNe la (SNe 2004dt, 2004ef, 2005M, and 2005cf) by 
using Hubble Space Telescope. One SN la in their sample, SN 
2004dt, displays a UV excess (the spectra reveal an excess in 
the 2900-3500 A wavelength range, compared with spectra of 
the other SN la events). In their study, the comparison object 
SN 2006X may also exhibit strong UV emission. The early 
UV emission may indicate the presence of a non-degenerate 
companion star in SN la progenitor systems. 

3.10. Early radio and X-ray emission of SNe la 

Circumstellar matter (CSM) provides a medium with which 
the SN ejecta can interact and produce radio synchrotron emis- 
sion. Many authors have searched for early radio emission from 
SNe la, but no detection has been made (Weiler et al., 1989; Eck 
et al., 1995, 2002). Hancock et al. (201 1) recently have used a 
stacking analysis of 46 archival Very Large Array observations 
by Panagia et al. (2006) to set upper limits on the radio emis- 
sion from SNe la in nearby galaxies. They gave an upper limit 
on the SN la peak radio luminosity of 1 .2 x 10 25 erg s _1 Hz T 1 at 



5 GHz, which implies an upper limit on the average companion 
stellar wind mass-loss rate of 1.3 x lO _7 M yr _1 before a SN 
explosion. Hancock et al. (201 1) argued that these limits chal- 
lenge expectations if the SN ejecta were encountering a CSM 
from the SD model. 

Aside from radio emission, the interaction of SN ejecta 
with the CSM can also produce X-ray emission. SN shock 
would run into CSM and heat it to high enough tempera- 
tures (~10 6 - 10 9 K), resulting in thermal X-rays (Chevalier, 
1990). Compared with radio emission, X-rays from SNe la re- 
sult from a different process and from different regions in the 
shocked CSM. Thus, it is a completely independent method 
to constrain progenitor model via detecting early X-ray emis- 
sion of SNe la. Russel and Immler (2012) recently consid- 
ered 53 SNe la observed by the Swift X-Ray Telescope. They 
gave an upper limit on the X-ray luminosity (0.2 - lOkeV) of 
1.7xl0 38 erg s _1 , which implies an upper limit on mass-loss rate 
of 1.1 x 10 _6 M o yr _1 x (v w )/(10kms _1 ), where v w is the wind 
speed for red supergiants that ranges from 5 to 25 km s _1 . Rus- 
sel and Immler (2012) claimed that these limits exclude mas- 
sive or evolved stars as the companions in progenitor systems 
of SNe la, but allow the possibility of MS and WD as the com- 
panion. 

According to the spin-up/spin-down model of SNe la sug- 
gested by Justham (2011) and Di Stefano et al. (2011), there 
is a delay between the time at which the WD has completed 
its mass-accretion and the time of the SN explosion. Since the 
matter ejected from the binary system during the mass-transfer 
has a chance to become diffuse, the SN explosion will occur in 
a medium with a density similar to that of typical regions of the 
interstellar medium. Therefore, the SD model may be compati- 
ble with the upper limits from SN la radio and X-ray detection. 

3.11. Pre-explosion images 

One of the methods to clarify SN la progenitor models is 
to directly detect the progenitor of a SN la in pre-explosion im- 
ages of the position where the SN occurred. Voss and Nelemans 
(2008) first studied the pre-explosion archival X-ray images at 
the position of the recent SN 2007on, and considered that its 
progenitor may be a WD + RG system. However, Roelofs et 
al. (2008) did not detect any X-ray source in images taken six 
weeks after SN 2007on's optical maximum and found an offset 
between the SN and the measured X-ray source position. Nele- 
mans et al. (2008) also obtained an ambiguous answer. Nielsen 
et al. (201 1) recently derived the upper limits of the X-ray lu- 
minosities from the locations of ten SNe la in nearby galaxies 
(<25 Mpc) before the explosions, most above a few 10 38 erg s" 1 
(for details see Fig. 1 of Nielsen et al., 201 1), which indicates 
that the progenitors of these SNe la were not bright supersoft 
X-ray sources shortly before they exploded as SNe la. How- 
ever, the upper limits are not constraining enough to rule out 
less bright supersoft X-ray progenitors (Nielsen et al., 2011). 
Future observations may shed light on the connection between 
SN la progenitors and X-ray emission. 

SN 201 lfe occurred in M101 at a distance of 6.4 Mpc is the 
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second closest SN la in the digital imaging era0 which was 
discovered by the Palomar Transient Factory survey less than a 
day after its explosion (Nugent et al., 2011a), and quickly fol- 
lowed up in many wavebands (Li et al., 2011; Nugent et al., 
201 lb; Smith et al., 201 1; Tammann and Reindl, 201 1; Patat et 
al., 201 lb; Liu et al., 2012; Horesh et al., 2012; Chomiuk et al., 
2012; Bloom et al., 2012; Brown et al., 2012a; Margutti et al., 
2012). Li et al. (201 1) used extensive historical imaging ob- 
tained at the location of SN 201 lfe to constrain the visible-light 
luminosity of the progenitor to be 10-100 times fainter than 
previous limits on other SN la progenitors. This result rules out 
luminous RG stars and most He stars as the mass donor star of 
this SN progenitor. These observations favour a scenario where 
the progenitor of SN 201 lfe accreted material either from WD, 
or via RLOF from a MS or subgiant companion. In a subse- 
quent work, Liu et al. (2012) also excluded its progenitor sys- 
tem with the most hottest photospheres by constraining X-ray 
properties prior to the SN explosion. 

Very recently, Horesh et al. (2012) set upper limits on both 
radio and X-ray emission from SN 201 lfe, excluding the pres- 
ence of a circumstellar matter from a giant donor star. Based on 
deep radio observations, Chomiuk et al. (2012) also excluded 
the presence of circumstellar matter. By using early optical and 
UV observations of SN 201 lfe, Nugent et al. (201 lb) excluded 
the presence of shock effects from SN ejecta hitting a compan- 
ion, and put a strict upper limit to the exploding star radius 
(<0.1 R ), thus providing a direct evidence that the progenitor 
is a compact star. A recent study by Bloom et al. (2012) also 
ruled out a MS star as the mass donor star and seem to favor 
a DD progenitor for SN 201 lfe (also see Brown et al., 2012a). 
We note that the spin-up/spin-down model potentially affects 
the conclusions above. 

3.12. Polarization ofSNe la 

Spectropolarimetry provides a direct probe of early time SN 
geometry, which is an important diagnostic tool for discriminat- 
ing among SN la progenitor systems and theories of SN explo- 
sion physics (see Livio and Pringle, 201 1). A hot young SN at- 
mosphere is dominated by the electron scattering that is highly 
polarizing. For an unresolved source with a spherical distribu- 
tion of scattering electrons, the directional components of the 
electric vectors of the scattered photons counteract exactly, re- 
sulting in zero net linear polarization. However, an incomplete 
cancelation will be derived from any asymmetry in the distribu- 
tion of the scattering electrons, or of absorbing material over- 
lying the electron-scattering atmosphere, which produces a net 
polarization (Leonard and Filippenko, 2005). 

SN asymmetry can therefore be measured via spectropo- 
larimetry, since asymmetric electron scattering leads to polar- 
ization vectors that do not cancel. Most normal SNe la are 
found to be spherically symmetric (a rather low polarization, 
<0.3%; Wang et al., 1996; Wang and Wheeler, 2008), but 
asymmetry has been detected at significant levels for a range 



The closest SN la in the digital imaging era is SN 1986G that exploded in 
NGC 5128 at a distance of ~4Mpc (Frogel et al., 1987). 



of SN la subclasses, e.g. sub-luminous SNe la with a contin- 
uum polarization about 0.3%-0.8% (Howell et al., 2001), and 
high-velocity (HV) SNe la with a high polarization about 2%, 
the spectra of which around maximum light are characterized 
by unusually broad and highly blueshifted absorption troughs 
in many line features (Leonard et al., 2005). Leonard et al. 
(2005) claimed that the following order emerges in terms of 
increasing strength of line-polarization features: normal/over- 
luminous SNe la < sub-luminous SNe la < HV SNe la. They 
argued that the most convincing explanation for the linear po- 
larization of all objects is partial obscuration of the photosphere 
by clumps of intermediate-mass elements forged in the SN ex- 
plosion. For a review of SN la polarimetric studies see Wang 
and Wheeler (2008). 

The explosion mechanism itself may produce asymmetry due 
to off-center explosion, and thus a polarization spectrum is ex- 
pected (Plewa et al., 2004; Kasen and Plewa, 2005). Thus, it is 
possible to obtain insight into the SN explosion physics with 
spectropolarimetry. Meanwhile, the progenitor systems may 
also cause the asymmetry. The SD model provides a natural 
way to produce the asymmetry. The existence of a compan- 
ion in the SD model may change the configuration of the SN 
ejecta (e.g. a cone-shaped hole shadowed by the companion), 
and thus a polarization spectrum is expected (Marietta et al., 
2000; Kasen et al., 2004; Meng and Yang, 2010b)Q In addi- 
tion, the DD model may also naturally result in an asymmetry 
of the distribution of SN ejecta. One relevant mechanism is the 
rapid rotation of a WD before a SN explosion, which leads to 
a change in the stellar shape. Another is that there may be a 
thick accretion disc around the CO WD, which may be an ori- 
gin of asymmetry in the configuration of the SN ejecta (e.g. 
Hillebrandt and Niemeyer, 2000). 

Livio and Pringle (201 1) argued that the nature of the corre- 
lation between the polarization and the observed SN la proper- 
ties can be used to distinguish between the SD and DD models. 
As a specific example, they considered possible correlations be- 
tween the polarization and the velocity gradient; a SN explosion 
is viewed from one pole it is seen as a high velocity gradient 
event at early phases with redshifts in late-time emission lines, 
while if it is viewed from the other pole it is seen as a low- 
velocity gradient event with blueshifts at late phases (Maeda et 
al, 2010). In the SD model, it is expected that the velocity gradi- 
ent is a two-valued function of polarization, with the largest and 
smallest values corresponding to essentially zero polarization. 
In the DD model, it is expected that the observed SN properties 
(i.e. velocity gradient) is a single-valued and monotonic func- 
tion of polarization. For details see Fig. 1 of Livio and Pringle 

(2011) . 

4. Related objects 

There are some objects that may be related to the progeni- 
tors and surviving companions of SNe la in observations, e.g. 

7 By using smoothed particle hydrodynamics simulations, Garcfa-Senz et al. 

(2012) studied the interaction of the hole, SN material and ambient medium. 
They concluded that the hole could remain open in the SNR for hundreds of 
years, suggesting the hole could affect its structure and evolution. 
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supersoft X-ray sources, cataclysmic variables, symbiotic sys- 
tems, single low-mass He WDs and hypervelocity He stars, etc. 

4.1. Supersoft X-ray sources 

Supersoft X-ray sources (SSSs) are one of the most promis- 
ing progenitor candidates of SNe la. Binaries in which steady 
nuclear burning takes place on the surface of the WDs have 
been identified with bright SSSs, discovered by the ROSAT 
satellite (van den Heuvel et al., 1992; Rappaport et al., 1994; 
Kahabka and van den Heuvel, 1997). Most of the known 
SSSs are located in the Large Magellanic Cloud, Small Mag- 
ellanic Cloud and M31. They typically emit 10 36 - 10 38 erg s -1 
in the form of very soft X-rays, peaking in the energy range 
20-100eV. 

van den Heuvel et al. (1992) proposed a model that the rel- 
atively massive WD sustains steady H-burning from a MS or 
subgiant donor star. They suggested that the mass-accretion 
occurs at an appropriate rate, in the range of 1.0 - 4.0 X 
lO~ 7 M yr _1 . Meanwhile, a WD + He star system has lumi- 
nosity around 10 37 - 10 38 ergs~' when the He-burning is sta- 
ble on the surface of the WD, which is consistent with that of 
observed from SSSs. Thus, WD + He star systems may also 
appear as SSSs before SN explosions (Iben and Tutukov, 1994; 
Yoon and Langer, 2003; Wang et al., 2009a). In addition, in 
the context of SSSs, the time that elapses between the double 
WD merger and the SN explosion is about 10 5 yr, and during 
this phase the merged object would look like as a SSS (with 
T ~ 0.5 - 1 x 10 6 K and L X -ra y ~ 10 37 ergs _1 ), which could 
provide a potential test for the DD model (Yoon et al., 2007; 
Voss and Nelemans, 2008). Note that the Galactic interstellar 
absorption and circumstellar matter may play an important role 
in the obscuration of X-rays. 

Recently, Di Stefano (2010a,b) called attention to the fact 
that in the galaxies of different morphological types there ex- 
ists a significant (up to 2 orders of magnitude) deficit of SSSs 
as compared with expectations based on SN la birthrates from 
the SD model. Gilfanov and Bogdan (2010) also obtained the 
same conclusion, based on the study of the luminosity of el- 
liptical galaxies in the supersoft X-ray range. However, these 
authors did not consider the binary evolution. A typical binary 
in the SD model undergoes three evolutionary stages in order of 
time before SN explosion, i.e. the wind phase, the supersoft X- 
ray source phase and the recurrent nova phase, since the mass- 
accretion rate decreases with time as the mass of the donor star 
decreases. The supersoft X-ray source phase is only a short 
time (e.g. a few hundred thousand years), since the SD progen- 
itor system spends a large part of lifetime in the wind phase or 
recurrent nova phase on its way to SN explosion (e.g. Han and 
Podsiadlowski, 2004; Meng et al., 2009; Wang et al., 2009a, 
2010a; Hachisu et al., 2010; Meng and Yang, 2011a). Lipunov 
et al. (2011) also considered that the theoretical SSS lifetimes 
and X-ray luminosities have been overestimated. 

4.2. Cataclysmic variables 

Cataclysmic variable stars (CVs) are stars that irregularly in- 
crease in brightness by a large factor, then drop back down to 



a quiescent phase (Warner, 1995). They consist of two com- 
ponent stars: a WD primary and a mass donor star. CVs are 
usually divided into several types, such as classical novae, re- 
current novae, nova-like variables, dwarf novae, magnetic CVs 
and AM CVns, etc (Warner, 1995). Among these subclasses of 
CVs, recurrent novae and dwarf novae are the most probable 
candidates of SN la progenitors. 

Recurrent novae have outbursts of about 4-9 magnitudes, 
and exhibit multiple outbursts at intervals of 10-80 years 
(Warner, 1995). They contain a massive WD and a relatively 
high mass-accretion rate (but below steady burning rate). The 
evolution of the outburst is very fast. Since the heavy element 
enhancement is not detected in recurrent novae, their WD mass 
is supposed to increase after each outburst. Additionally, nova 
outbursts require a relatively high mass-accretion rate onto a 
massive WD to explain the recurring nova outbursts. Thus, 
these objects become some of the most likely candidates of 
SN la progenitors (Starrfield et al., 1985; Hachisu and Kato, 
2001). However, this class of objects are rare, with ten Galactic 
recurrent novae, two in the Large Magellanic Cloud and a few 
in M31. Recurrent novae and SSSs differ in the mass-accretion 
rate from a mass donor star onto the WD; SSSs have steady nu- 
clear burning on the surface of the WD, while recurrent novae 
happen at rates that allow shell flashes. 

By modeling the decline of the outburst light curves of some 
recurrent novae (T CrB, RS Oph, V745 Sco and V3890 Sgr), 
Hachisu and Kato (2001) suggested that these WDs are ap- 
proaching the Ch mass and will produce SNe la. Recurrent 
nova systems like RS Oph have been proposed as possible SN la 
progenitors, based on the high mass of the accreting WD. Patat 
et al. (201 la) investigated the circumstellar environment of RS 
Oph and its structure, suggesting that the recurrent eruptions 
might create complex structures within the material lost by the 
donor star. This may establish a strong link between RS Oph 
and the progenitor system of SN 2006X, for which similar fea- 
tures have been detected. 

Recurrent nova U Sco contains a WD of M W d = 1.55 ± 
0.24 M and a secondary star with M 2 = 0.88 + 0.17 M or- 
biting with a period P or b ~ 0.163 day (Thoroughgood et al., 
2001). The high mass of the WD implies that U Sco is a strong 
progenitor candidate of a SN la (Thoroughgood et al., 2001; 
also see Hachisu et al., 2000). However, the nebular spectra 
of U Sco displays that the relative abundance of [Ne/O] is 1 .69, 
which is higher than that of the typical [Ne/O] abundance found 
in classical novae from CO WDs and suggests that U Sco has a 
O-Ne-Mg WD (Mason, 201 1). Thus, U Sco may not explode as 
a SN la but rather collapse to a neutron star by electron capture 
on 24 Mg. 

Dwarf novae have multiple outbursts ranging in brightness 
from 2 to 5 magnitudes, and exhibit intervals from days to 
decades. The lifetime of an outburst is typically from 2 to 20 
days and is related to the outburst interval. Dwarf nova out- 
bursts are usually attributed to the release of gravitational en- 
ergy resulted from an instability in the accretion disk or by sud- 
den mass-transfer via the disk (Warner, 1995). Observationally, 
there are a number of dwarf novae in which the WD is about 
1 M (e.g. GK Per, SS Aur, HL CMa, U Gem, Z Cam, SY Cnc, 
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OY Car, TW Vir, AM Her, SS Cyg, RU Peg, GD 552 and IP 
Peg, etc). The secondaries of these WD binaries are K or M 
stars (<1 M Q ). A few of these systems with early K type sec- 
ondaries may have the WD mass close to the Ch mass. It has 
been suggested that the mass-accretion rate onto a WD during 
a dwarf nova outbursts can be sufficiently high to allow steady 
nuclear burning of the accreted matter and growth of the WD 
mass (King et al., 2003; Xu and Li, 2009; Wang et al., 2010a; 
Meng and Yang, 2010a). However, whether dwarf nova out- 
bursts can increase the mass of a WD close to Ch mass is still a 
problem (e.g. Hachisu et al., 2010). 

4.3. Symbiotic systems 

Symbiotic systems are long-period binaries, consisting of a 
RG and a hot object that is usually a WD (Truran and Cameron, 
1971). The hot object accretes and burns material from the RG 
star via stellar wind in most cases, but could also be RLOF in 
some cases. They usually show strong emission lines from sur- 
rounding circumstellar material ionized by the hot component, 
and low temperature absorption features from the RG. Symbi- 
otic systems are essential to understand the evolution and inter- 
action of detached and semi-detached binaries. There are two 
distinct subclasses of symbiotic stars, i.e. the S-type (stellar) 
with normal RG stars and orbital periods of about 1-15 years, 
and the D-type (dusty) with Mira primaries usually surrounded 
by a warm dust shell and orbital periods longer than 10 years. 
Symbiotic stars are thus interacting binaries with the longest 
orbital periods. Tang et al. (2012) recently found a peculiar 
symbiotic system J0757 that consists of an accreting WD and a 
RG. In quiescent phase, however, it doesn't show any signature 
of "symbiotic". Thus, it is a missing population among sym- 
biotic systems, which may contribute to a significant fraction 
of SN la. Moreover, this object showed a 10 year flare in the 
1940s, possibly from H-shell burning on the surface of the WD 
and without significant mass-loss. Therefore, the WD could 
grow effectively. 

The presence of both the accreting WD and the RG star 
makes symbiotic binaries a promising nursery for the produc- 
tion of SNe la. However, due to the low efficiency of matter ac- 
cumulation by a WD accreting material from the stellar wind, 
SN la birthrate from these symbiotic systems is relative low 
(e.g. Yungelson and Livio, 1998). 

4.4. Single low-mass He WDs 

The existence of a population of single low-mass He WDs 
(LMWDs; <0.45 M©) is supported by some recent observations 
(e.g. Marsh et al., 1995; Kilic et al., 2007). However, it is still 
unclear how to form single LMWDs. It has been suggested that 
single LMWDs could be produced by single old metal-rich stars 
that experience significant mass-loss before the central He flash 
(Kalirai et al., 2007; Kilic et al., 2007). However, the study of 
the initial-final mass relation for stars by Han et al. (1994) im- 
plied that only LMWDs with masses larger than 0.4 M© might 
be produced from such a single star scenario, even at high 
metallicity environment (Meng et al., 2008). Thus, it would 
be difficult to conclude that single stars can produce LMWDs 
of -0.2 M Q . 



Justham et al. (2009) inferred an attractive formation sce- 
nario for single LMWDs, which could be formed in binaries 
where their companions have exploded as SNe la. Wang and 
Han (2010d) recently found that the surviving companions of 
the old SNe la from the WD + MS and WD + RG channels 
have low masses, providing a possible way to explain the for- 
mation of the population of single LMWDs (see also Meng and 
yang, 2010c). Conversely, the observed single LMWDs may 
provide evidence that at least some SN la explosions have oc- 
curred with non-degenerate donors (such as MS or RG donors). 
We note that Nelemans and Tauris (1998) also proposed an al- 
ternative scenario to form single LMWDs from a solar-like star 
accompanied by a massive planet, or a brown dwarf, in a rela- 
tively close binary orbit. 

4.5. Hypervelocity stars 

In recent years, hypervelocity stars (HVSs) have been ob- 
served in the halo of the Galaxy. HVSs are stars with velocities 
so high that they are able to escape the gravitational pull of the 
Galaxy. However, it is still not clear how to form HVSs (for 
a review see Tutukov and Fedorova, 2009). It has been sug- 
gested that such HVSs can be formed by the tidal disruption 
of a binary through interaction with the super-massive black 
hole (SMBH) at the Galactic center (GC) (Hills, 1988; Yu and 
Tremaine, 2003; Zhang et al., 2010). 

The first three HVSs have only recently been discovered 
serendipitously (e.g. Brown et al., 2005; Hirsch et al., 2005; 
Edelmann et al., 2005). Up to now, about 17 confirmed HVSs 
have been discovered in the Galaxy (Brown et al., 2009; Tillich 
et al., 2009), most of which are B-type stars, probably with 
masses ranging from 3 to 5 M© (Brown et al., 2005, 2009; Edel- 
mann et al., 2005). The HVS B-type stars are demonstrated 
short-lived B-type stars at 50-100kpc distances that are signif- 
icantly unbound based on radial velocity alone. Their observed 
properties (ages, flight times, latitude distribution) are consis- 
tent with the Galactic center ejection scenario (Brown et al., 
2012b). One HVS, HE 0437-5439, is known to be an apparently 
normal early B-type star. Edelmann et al. (2005) suggested that 
the star could have originated in the Large Magellanic Cloud, 
since it is much closer to this galaxy (~18kpc) than to the GC 
(see also Przybilla et al., 2008). Li et al. (2012) recently re- 
ported 13 metal -poor F-type HVS candidates which are selected 
from 370,000 stars of the data release 7 of the Sloan Digital 
Sky Survey. With a detailed analysis of the kinematics of these 
stars, they claimed that seven of them were likely ejected from 
the GC or the Galactic disk, four neither originated from the GC 
nor the Galactic disk, and the other two were possibly ejected 
from either the Galactic disk or other regions. 

At present, only one HVS, US 708, is an extremely He-rich 
sdO star in the Galactic halo, with a heliocentric radial velocity 
of +708 ± 15 km/s. Hirsch et al. (2005) speculated that US 708 
was formed by the merger of two He WDs in a close binary 
induced by the interaction with the SMBH in the GC and then 
escaped. Recently, Perets (2009) suggested that US 708 may 
have been ejected as a binary from a triple disruption by the 
SMBH, which later on evolved and merged to form a sdO star. 
However, the evolutionary lifetime of US 708 is not enough if 
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it originated from the GC. Wang and Han (2009) found that the 
surviving companions from the He star donor channel have a 
high spatial velocity (>400km/s) after a SN explosion, which 
could be an alternative origin for HVSs, especially for HVSs 
such as US 708 (see also Justham et al., 2009). Considering 
the local velocity nearby the Sun (~220km/s), Wang and Han 
(2009) found that about 30% of the surviving companions may 
be observed to have velocity above 700km/s. In addition, a 
SN asymmetric explosion may also enhance the velocity of the 
surviving companion. Thus, a surviving companion star in the 
He star donor channel may have a high velocity like US 708. 

5. Origin of SN la diversity 

SNe la have been successfully used as cosmological distance 
candles, but there exists spectroscopic diversity among SNe la 
that is presently not well understood, nor how this diversity is 
linked to the properties of their progenitors (e.g. Branch et al., 
1995; Livio, 2000). When SNe la are applied as distance in- 
dicators, the Phillips relation is adopted (i.e. the luminosity- 
width relation; brighter SNe la have wider light curves), which 
implies that SN la luminosity is mainly determined by one pa- 
rameter. In an attempt to quantify the rate of spectroscopi- 
cally peculiar SNe la in the existing observed sample, Branch 
et al. (1993) compiled a set of 84 SNe la and found that about 
83% - 89% of the sample are normal. According to the study 
of Li et al. (2001), however, only 64% + 12% of the observed 
SNe la are normal in a volume-limited searchJI The total rate 
of peculiar SNe la could be as high as 36% + 9%; the rates are 
16% + 7% and 20% ± 7% for SN 1991bg-like objects and SN 
1991T-fike objects, respectively. SN 1991bg-like objects both 
rise to their maximum and decline more quickly, and are sub- 
luminous relative to normal SNe la, whereas SN 1991T-like ob- 
jects both rise to their maximum and decline more slowly, and 
are more luminous relative to normal SNe la. These two types 
of peculiar events obey the luminosity-width relation. However, 
a subset of SNe la apparently deviate from the luminosity-width 
relation, e.g. some were observed with exceptionally high lumi- 
nosity or extremely low luminosity, which may have progeni- 
tors with masses exceeding or below the standard Ch mass limit 
(e.g. Howell et al., 2006; Foley et al., 2009). This implies that 
at least some SNe la can be produced by a variety of different 
progenitor systems, and probably suggests that SN la luminos- 
ity is not the single parameter of the light curve shape. 

It has been suggested that the amount of 56 Ni formed during 
a SN la explosion dominates its maximum luminosity (Arnett, 
1982), but the origin of the variation of the amount of 56 Ni for 



There is increasing evidence showing that even the normal SNe la exhibit 
diversity in their spectral features (e.g. Branch et al., 2009; Wang et al., 2009c; 
Blondin et al., 2012). Wang et al. (2009c) investigated 158 relatively normal 
SNe la by dividing them into two groups in terms of the expansion velocity 
inferred from the absorption minimum of the Sill /16355 line around maximum 
light. They claimed that, one group "Normal" consists of SNe la with an av- 
erage expansion velocity 10, 600 ± 400km/s, but another group "HV" consists 
of objects with higher velocities ~1 1, 800km/s. The HV SNe la are found to 
prefer a smaller extinction ratio Ry (relative to the Normal ones), which might 
suggest the presence of circumstellar material (see Sect. 3.7). 



different SNe la is still unclear (the derived Ni masses for dif- 
ferent SNe la could vary by a factor of ten; Wang et al., 2008a). 
Many efforts have been paid to solve this problem. Umeda et 
al. (1999) suggested that the average ratio of carbon to oxy- 
gen (C/O) of a WD at the moment of a SN explosion is the 
dominant parameter for the Phillips relation, i.e. the higher the 
C/O ratio, the larger the amount of 56 Ni, and then the higher the 
maximum luminosity (see also Meng and Yang, 201 lb). How- 
ever, 3D simulations by Ropke and Hillebrandt (2004) suggest 
that different C/O ratios have a negligible effect on the amount 
of 56 Ni produced. At present, the studies from the explosion 
models of SNe la indicate that the number of ignition points at 
the center of WDs or the transition density from deflagration to 
detonation dominates the production of 56 Ni, and consequently 
the maximum luminosity (e.g. Hillebrandt and Niemeyer, 2000; 
Hoflich et al., 2010; Kasen et al., 2010). 

It was claimed that the ignition intensity (the number of igni- 
tion points) in the center of WDs is a useful parameter in inter- 
preting the Phillips relation (Hillebrandt and Niemeyer, 2000). 
Based on the SD model, Lesaffre et al. (2006) carried out a sys- 
tematic study of the sensitivity of carbon ignition conditions for 
the Ch mass WDs on various properties, and claimed that the 
central density of a WD at the carbon ignition may be the ori- 
gin of the scatter of the maximum luminosity. This suggestion 
was further supported by detailed multi-dimensional numerical 
simulations of SN explosions (Rrueger et al., 2010). We note 
that the WD cooling time before mass-accretion is less than 
1 Gyr in the simulations of Lesaffre et al. (2006) and Krueger 
et al. (2010). However, there are SNe la with the delay times 
~10 Gyr in observations. The WDs with such a long cooling 
time may become more degenerate before the onset of the mass- 
accretion phase. Some other processes, such as carbon and oxy- 
gen separation or crystallization, may occur and dominate the 
properties of the CO WD (Fontaine et al., 2001). How the ex- 
tremely degenerate conditions affect the properties of SNe la 
still remains unclear. The suggestion of Lesaffre et al. (2006) 
should be checked carefully under extremely degenerate con- 
ditions. Adopting the WD mass-accretion process in Lesaffre 
et al. (2006), Chen et al. (2012) recently studied the evolution 
of various CO WDs from the onset of mass-accretion to carbon 
ignition at Ch mass limit. The study shows that the carbon ig- 
nition generally occurs at the center for hot low-mass CO WDs 
but off-center for cool massive ones, which may provide more 
information for the explosion models of SNe la. 

Some numerical and synthetical results showed that the 
metallicity may have an effect on the final amount of 56 Ni, and 
thus the maximum luminosity of SNe la (Timmes et al., 2003; 
Podsiadlowski et al., 2006; Bravo et al., 2010). There is also 
some other evidence of the correlation between the properties 
of SNe la and metallicity from observations (e.g. Branch and 
Bergh, 1993; Hamuy et al., 1996; Wang et al, 1997; Gallagher 
et al., 2008; Sullivan, 2006; Howell et al., 2009a; Sullivan et 
al., 2010). Podsiadlowski et al. (2006) introduced metallicity 
as a second parameter that affects the light curve shape. For 
a reasonable range of metallicity, this may account for the ob- 
served spread in the Phillips relation. Since metallicity in the 
Universe has evolved with time, this introduces an undesirable 
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evolutionary effect in the SN la distance method, which could 
mimic the effect of an accelerating Universe. We also note that 
Maeda et al. (2010) argued that the origin of spectral evolu- 
tion diversity in SNe la can be understood by an asymmetry in 
the SN explosion combined with the observer's viewing angle. 
Moreover, Parrent et al. (201 1) investigated the presence of CII 
/16580 in the optical spectra of 19 SNe la. Most of the objects 
in their sample that exhibit CII /16580 absorption features are 
of the low-velocity gradient subtype. This study indicates that 
the morphology of carbon-rich regions is consistent with either 
a spherical distribution or a hemispheric asymmetry, support- 
ing the idea that SN la diversity may be a result of off-center 
ignition coupled with observer's viewing angle. 

6. Impacts of SN la progenitors on some fields 

The identification of SN la progenitors also has important 
impacts on some other astrophysical fields, e.g. cosmology, the 
evolution of galaxies, SN explosion models and binary evolu- 
tion theories, etc (e.g. Branch et al., 1995; Livio, 2000). 

Cosmology. It is feasible to improve SNe la as mature 
cosmological probes, since the dominant systematic errors are 
clear, which include photometric calibration, selection effects, 
reddening and population-dependent differences, etc. In the 
next decade, SNe la are proposed to be cosmological probes 
for testing the evolution of the dark energy equation of state 
with time (Howell et al., 2009b). The use of SNe la as one 
of the main ways to determine the Hubble constant (Hq) and 
cosmological parameters (e.g. £1 M and Q A ; Riess et al., 1998; 
Perlmutter et al., 1999), requires our understanding of the evo- 
lution of the luminosities and birthrates of SNe la with cosmic 
epoch. Both of these depend on the nature of their progeni- 
tors. Meanwhile, the evolution of the progenitor systems or a 
changing mix of different progenitors may bias cosmological 
inferences. For a recent review of this field see Howell (201 1). 

Galaxy evolution. Aside from cosmology, the evolution of 
galaxies depends on the radiative, kinetic energy, nucleosyn- 
thetic outputs (e.g. Kauffmann et al., 1993; Liu et al., 2012a) 
and the birthrates of SNe la with time, which all depend on the 
nature of the progenitor systems. SNe la are also laboratories 
for some extreme physics, e.g. they are accelerators of cos- 
mic rays and as sources of kinetic energy in galaxy evolution 
processes (e.g. Helder et al., 2009; Powell et al., 2011). Es- 
pecially, SNe la regulate galactic and cluster's chemical evolu- 
tion. Due to the main contribution of iron to their host galaxies, 
SNe la are a key part of our understanding of galactic chemi- 
cal evolution (e.g. Greggio and Renzini, 1983; Matteucci and 
Greggio, 1986). The existence of young and old populations of 
SNe la suggested by recent observations may have an impor- 
tant effect on models of galactic chemical evolution, since they 
would return large amounts of iron to the interstellar medium 
either much earlier or much later than previously thought. 

Explosion models. SNe la provide natural laboratories for 
studying the physics of hydrodynamic and nuclear processes 
with extreme conditions. The link between the progenitor mod- 
els and the explosion models is presently one of the weak- 
est points in our understanding of SNe la (Hillebrandt and 



Niemeyer, 2000). Due to some uncertainties that still exist in 
the SN explosion mechanism itself, a knowledge of the initial 
conditions and the distribution of matter in the environment of 
the exploding star is essential for our understanding of SN ex- 
plosion, e.g. the ignition density may depend on the initial WD 
mass, the age of the progenitor, the metallicity and the treat- 
ment of rotation in the progenitor. Moreover, different pro- 
genitor models may lead to different WD structures before SN 
explosion. Lu et al. (2011) recently studied the properties of 
the Tycho's SNR. They estimated the parameters of the binary 
system before the SN explosion, which may shed lights on the 
possible explosion models. 

Binary evolution theories. The identification of SN la 
progenitors, coupled with observationally determined SN la 
birthrates and delay times will help to place meaningful con- 
straints on some theories of binary evolution, e.g. the mass- 
transfer between two stars, the mass-accretion efficiency of 
WDs, etc (e.g. Hachisu et al., 1996; Han and Podsiadlowski, 
2004; Wang et al., 2009a). Especially, it is possible that the 
CE efficiency parameter may be constrained (e.g. Meng et al, 
201 1), which is important in binary evolution and BPS studies. 

7. Summary 

In this article, various progenitor models proposed in the lit- 
eratures are reviewed, including some variants of SD and DD 
models. We addressed some observational ways to test the cur- 
rent progenitor models and introduced some observed objects 
that may be related to the progenitors and the surviving com- 
panion stars of SNe la. We also discussed the impacts of SN la 
progenitors on some fields. The origin of the observed SN la 
diversity is still unclear. It seems likely that SNe la can be pro- 
duced by a variety of different progenitor systems, perhaps ex- 
plaining part of the observed diversity. SN asymmetric explo- 
sion coupled with observer's viewing angle may also produce 
the diversity. Additionally, the metallicity of progenitors may 
be a second parameter that affects the light curve shape of SNe 
la. 

At present, the SD model is the most widely accepted SN la 
progenitor model. The advantages of this model can be sum- 
marized as follows: 

(1) The SD model is in excellent agreement with the observed 
light curves and spectroscopy of SNe la, and this model 
may explain the similarities of most SNe la. 

(2) Observationally, there is increasing evidence indicating that 
some SNe la may come from the SD model (e.g. the signa- 
tures of gas outflows from some SN la progenitor systems, 
the wind-blown cavity in SN remnant, and the early optical 
and UV emission of SNe la, etc). In addition, the SD model 
may be compatible with some recent observations (e.g. the 
lack of H or He seen in nebular spectra of SNe la, and the 
upper limits from SN la radio and X-ray detection, etc) by 
considering the spin-down time. 

(3) There are some SD progenitor candidates in observations, 
e.g. supersoft X-ray sources, recurrent novae, dwarf novae 
and symbiotic systems, etc. Meanwhile, a number of high 
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mass WDs that have been accreting from a non-degenerate 
companion star have been found. 

(4) The observed single low-mass He WDs and hypervelocity 
He stars may be explained by the surviving companion stars 
predicted in the SD model. 

(5) SNe la with long delay times can be understood by the WD 
+ MS and WD + RG channels. In contrast, SNe la with 
short delay times may consist of systems with a He donor 
star in the WD + He channel, or even a massive MS donor 
star in the WD + MS channel. 

(6) Besides the DD model, these observed super-luminous 
SNe la can also be produced by the SD model by consider- 
ing the effects of rapid differential rotation on the accreting 
WD. 

However, the SD model is still suffering some problems from 
both theoretically and observationally that need to be resolved: 

(1) The optically thick wind assumption, widely adopted in the 
studies of the SD model, is in doubt for very low metallic - 
ity; the low-metallicity threshold for SNe la predicted by 
theories has not been found in observations. 

(2) It is still difficult to reproduce the observed birth rates and 
delay times of SNe la. This suggests that we need a better 
understanding of mass-accretion onto WDs. 

(3) There is still no conclusive proof that any individual object 
is the surviving companion star of a SN la, which is pre- 
dicted by the SD model. A likely surviving companion star 
for the progenitor of Tycho's SN has been identified, but 
the claim is still controversial. 

Although a DD merger is thought to experience an accretion- 
induced collapse rather than a thermonuclear explosion, any 
definitive conclusion about the DD model is currently prema- 
ture: 

(1) There are some parameter ranges in which the accretion- 
induced collapse can be avoided. Recent simulations in- 
dicate that the violent mergers of two massive WDs can 
closely resemble normal SN la explosion with the assump- 
tion of the detonation formation as an artificial parameter, 
although these mergers may only contribute a small frac- 
tion to the observed population of normal SNe la. 

(2) This model can naturally reproduce the observed birthrates 
and delay times of SNe la and may explain the formation 
of some observed super-luminous SNe la. 

(3) This model can explain the lack of H or He seen in the 
nebular spectra of SNe la. 

(4) Recent observational studies of SN 201 lfe seem to favor a 
DD progenitor. In addition, there is no signal of a surviving 
companion star from the central region of SNR 0509-67.5 
(the site of a SN la explosion whose light swept Earth about 
400 years ago), which may indicate that the progenitor for 
this particular SN la is a DD system. 

(5) Some observed double WD systems may have the total 
mass larger than the Ch mass, and possibly merge within 
the Hubble-time, although there are not enough double WD 
systems to reproduce the observed SN la birthrates in the 
context of the DD model. 



Some variants of the SD and DD models have been proposed 
to explain the observed diversity of SNe la: 

(1) The sub-luminous 1991bg-like objects may be explained 
by the sub-Ch mass model. 

(2) The unusual properties of 2002ic-like objects can be under- 
stood by the delayed dynamical instability model. 

(3) The spin-up/spin-down model may provide a route to ex- 
plain the similarities and the diversity observed in SNe la. 

(4) The core-degenerate model could form a massive WD with 
super-Ch mass that might explode as a super-luminous 
SN la. 

(5) The collisions between two WDs in dense environments 
could also potentially lead to sub-luminous SN la explo- 
sions. 

(6) The mechanism of WDs exploding near black holes is also 
a potential progenitor model for thermonuclear runaway, 
despite of the expected low rate when a WD passes near a 
black hole. 

To set further constraints on SN la progenitor models, large 
samples of SNe la with well-observed light curves and spec- 
troscopy in nearby galaxies are required to establish the connec- 
tion of SN la properties with the stellar environments of their 
host galaxies. Many new surveys from ground and space have 
been proposed to make strides in SN la studies , e.g. Palomar 
Transient Factory, Skymapper, La Silla QUEST, Pan-STARRS, 
the Dark Energy Survey, Large Synoptic Survey Telescope, the 
Joint Dark Energy Mission and the Gaia Astrometric Mission, 
etc (Howell et al., 2009; Altavilla et al., 2012). These surveys 
will allow comparisons via large SN la subsamples, and start 
to connect SN la progenitors with the observed features of SN 
explosions themselves, and thus to unveil the nature of SN la 
progenitors. 
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